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Abstract 


We solve for the cosmological perturbations in a five-dimensional background consisting of two 
separating or colliding boundary branes, as an expansion in the collision speed V divided by the 
speed of light c. Our solution permits a detailed check of the validity of four-dimensional effective 
theory in the vicinity of the event corresponding to the big crunch/big bang singularity. We show 
that the four-dimensional description fails at the first nontrivial order in ( V/c ) 2 . At this order, 
there is nontrivial mixing of the two relevant four-dimensional perturbation modes (the growing 
and decaying modes) as the boundary branes move from the narrowly-separated limit described by 
Kaluza-Klein theory to the well-separated limit where gravity is confined to the positive-tension 
brane. We comment on the cosmological significance of the result and compute other quantities of 
interest in five-dimensional cosmological scenarios. 
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I. INTRODUCTION 


One of the most striking implications of string theory and M theory is that there are extra 
spatial dimensions whose size and shape determine the particle spectrum and couplings 
of the low energy world. If the extra dimensions are compact and of fixed size R, their 
existence results in a tower of Kaluza-Klein massive modes whose mass scale is set by Hr 1 . 
Unfortunately, this prediction is hard to test if the only energy scales accessible to experiment 
are much lower than R^ 1 . At low energies, the massive modes decouple from the low energy 
effective theory and are, for all practical purposes, invisible. Therefore, we have no means 
of checking whether the four-dimensional effective theory observed to fit particle physics 
experiments is actually the outcome of a simpler higher-dimensional theory. 

The one situation where the extra dimensions seem bound to reveal themselves is in 
cosmology. At the big bang, the four-dimensional effective theory (Einstein gravity or its 
stringy generalization) breaks down, indicating that it must be replaced by an improved 
description. There are already suggestions of improved behavior in higher-dimensional string 
and M theory. If matter is localized on two boundary branes, the matter density remains 
finite at a brane collision even though this moment is, from the perspective of the four¬ 
dimensional effective theory, the big bang singularity [1, y, dj. Likewise, the equations of 
motion for fundamental strings are actually regular at t — 0 in string theory, in the relevant 
background solutions (J, 0] . 

In this paper, we shall not study the singularity itself. Instead, we will study the behavior 
of higher-dimensional gravity as the universe emerges from a brane collision. Our particular 
concern is to determine the extent to which the four-dimensional effective theory accurately 
captures the higher-dimensional dynamics near the big bang singularity. 

The model we study is the simplest possible model of braneworld gravity. It consists of 
two empty Z 2 -branes (or orbifold planes) of opposite tension separated by a five-dimensional 


bu 

L 


k with negative cosmological constant, corresponding to an anti-de Sitter (AdS) radius 
. Many works have been devoted to obtaining exact or approximate solutions for this 


model, for static or moving branes 


ed to obtaining e 

a 0 AH Oil- 


Our methods have much in common 


with these earlier works, in particular the idea that, when the branes move slowly, the four¬ 
dimensional effective theory works well. However, our focus and goals are rather different. 
When the two boundary branes are very close to one another, the warping of the hve- 
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dimensional bulk and the tension of the branes become irrelevant. In this situation, the 
low energy modes of the system are well-described by a simple Kaluza-Klein reduction from 
five to four dimensions, i.e., gravity plus a scalar held (the Z 2 projections eliminate the 
gauge held zero mode). We shall verify this expectation. However, when the two branes 
are widely separated, the physics is quite different. In this regime, the warping of the 
bulk plays a key role, causing the low energy gravitational modes to be localized on the 


positive-tension brane yjj. Ilf 


14 j . The four-dimensional effective theory describing this 


new situation is nevertheless identical, consisting of Einstein gravity and a scalar held, the 
radion, describing the separation of the two branes. 

In this paper, we study the transition between these two regimes - from the naive Kaluza- 
Klein reduction to localized Randall-Sundrum gravity - at finite brane speed. In the two 
asymptotic regimes - the narrowly-separated brane limit and the widely-separated limit - 
the cosmological perturbation modes show precisely the behavior predicted by the four¬ 
dimensional effective theory. There are two massless scalar perturbation modes; in longitu¬ 
dinal gauge, and in the long wavelength [k —> 0) limit, one mode is constant and the other 
decays as , where 1 4 is the conformal time. In the four-dimensional description, these two 
perturbation modes are entirely distinct: one is the curvature perturbation mode; the other 
is a local time delay to the big bang. However, we shall show that in the five-dimensional 
theory, at first nontrivial order in the speed of the brane collision, the two modes mix. If, for 
example, one starts out in the t im e delay mode at small £4, one ends up in a mixture of the 
time delay and curvature perturbation modes as t± —> 00. Thus the two cosmological pertur¬ 
bation modes - the growing and decaying adiabatic modes - mix in the higher-dimensional 
braneworld setup, a phenomenon which is prohibited in four dimensions. 

The mode-mixing occurs as a result of the qualitative change in the nature of the low 
energy modes of the system. At small brane separations the low energy modes are nearly 
uniform across the extra dimension. Yet as the brane separation becomes larger than the 
bulk warping scale, the low energy modes become exponentially localized on the positive- 
tension brane. If the branes separate at finite speed, the localization process fails to keep pace 
with the brane separation and the low energy modes do not evolve adiabatically. Instead, 
they evolve into a mixture involving higher Kaluza-Klein modes, and the four-dimensional 
effective description fails. 

As mentioned, the mixing we see between the two scalar perturbation modes would be 
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prohibited in any local four-dimensional effective theory consisting of Einstein gravity and 
matter fields, no matter what the matter fields were. Therefore the mixing is a truly five¬ 
dimensional phenomenon, which cannot be modeled with a local four-dimensional effective 
theory. There is an independent argument against the existence of any local four-dimensional 
description of these phenomena. In standard Kaluza-Klein theory, it is well known that the 


entire spectrum of massive modes is actually spin two Q|. Yet, despite many attempts, no 
satisfactory Lagrangian description of massive, purely spin two fields has ever been found 
Ifl Il7 |. Again, this suggests that one should not expect to describe the excitation of the 
higher Kaluza-Klein modes in terms of an improved, local, four-dimensional effective theory. 

The system we study consists of two branes emerging from a collision. In this situation, 
there are important simplifications which allow us to specify initial data rather precisely. 
When the brane separation is small, the fluctuation modes neatly separate into light Kaluza- 
Klein zero modes, which are constant along the extra dimension, and massive modes with 
nontrivial extra-dimensional dependence. Furthermore, the brane tensions and the bulk 
cosmological constant become irrelevant at short distances. It is thus natural to specify initial 
data which map precisely onto four-dimensional fields in the naive dimensionally-reduced 
theory describing the limit of narrowly-separated branes. With initial data specified this 
way, there are no ambiguities in the system. The two branes provide boundary conditions 
for all time and the five-dimensional Einstein equations yield a unique solution, for arbitrary 
four-dimensional initial data. 

Our main motivation is the study of cosmologies in which the big bang was a brane col¬ 


lision, such as the cyclic model [ 1 ], Here, a period of dark energy domination, followed by 
slow contraction of the fifth dimension, renders the branes locally flat and parallel at the 
collision. During the slow contraction phase, growing, adiabatic, scale-invariant perturba¬ 
tions are imprinted on the branes prior to the collision. However, if the system is accurately 
described by four-dimensional effective theory throughout, then, as a number of authors 
have noted U Q Q 0 , 122I 23], there is an apparent roadblock to the passage of the 
scale-invariant perturbations across the bounce. Namely, it is hard to see how the growing 
mode in the contracting phase, usually described as a local time delay, could match onto 
the growing mode in the expanding phase, usually described as a curvature perturbation. In 
this paper, we show that the four-dimensional effective theory fails at order ( Vjc ) 2 , where 
V is the collision speed and c is the speed of light. The four-dimensional description works 
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well when the branes are close together, or far apart. However, as the branes move from one 
regime to the other, the two four-dimensional modes mix in a nontrivial manner. 

The mixing we find is an explicit demonstration of a new physical principle: namely, 
the approach of two boundary branes along a fifth dimension produces physical effects that 
cannot properly be modeled by a local four-dimensional effective theory. In this paper, we 
deal with the simplest case involving two empty boundary branes separated by a bulk with 
a negative cosmological constant. For the cyclic model, the details are more complicated 


24]. The bulk possesses an additional bulk stress AT|, associated with the inter-brane 


force, that plays a vital role in converting a growing mode corresponding to a pure time 
delay perturbation into a mixture of time delay and curvature modes on the brane. We also 
explain how subsequent five-dimensional effects cause cause the curvature on the brane to 
feed into the adiabatic growing mode perturbation in the four-dimensional effective theory 
after the bang. The details are different, but the general principle is the same as in the case 
considered in this paper. One must go beyond the four-dimensional effective theory and 
consider the full five-dimensional theory to compute properly the evolution of perturbations 
before and after a brane collision. 

Even though our main concern is with cyclic/ekpyrotic models, our methods are likely to 
be more broadly applicable in braneworld models. (See e 


■g. 0 


26 


27] for reviews). Our 


methods may be extended, for example, to models with better motivation from fundamental 
theory, such as heterotic M theory |2Sl 123] • One may also include matter on the branes. 
Another interesting application would be to study the evolution of a four-dimensional black 
hole in the limit of narrowly-separated branes, i.e., a black string, as the two branes separate 
and the Gregory-Laflamme instability appears. 

The outline of this paper is as follows. In Section II we provide an overview of our three 
solution methods. In Section III we solve for the background and cosmological perturbations 
using a series expansion in time about the collision. In Section IV we present an improved 
method in which the dependence on the fifth dimension is approximated using a set of higher- 
order Dirichlet or Neumann polynomials. In Section V we develop an expansion about the 
small-(V/c) scaling solution, before comparing our results with those of the four-dimensional 
effective theory in Section VI. We conclude with a discussion of mode-mixing in Section VII. 
Detailed explicit solutions may be found in the Appendices, and the Mathematica code 
implementing our calculations is available online 
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II. THREE SOLUTION METHODS 


In this section, we review the three solution methods employed, noting their comparative 
merits. For the model considered here, with no dynamical bulk helds, there is a Birkhoff- 
like theorem guaranteeing the existence of coordinates in which the bulk is static. It is 
easy to solve for the background in these coordinates. However, the motion of the branes 
complicates the Israel matching conditions rendering the treatment of perturbations difficult. 
It is preferable to choose a coordinate system in which the branes arc located at fixed spatial 
coordinate y = ±?/o> and the bulk evolves with time. 

We shall employ a coordinate system in which the five-dimensional line element for the 
background takes the form 

ds 2 = n 2 (t, y){—dt 2 + t 2 dy 2 ) + & 2 (t, y) da; 2 , (1) 

where y parameterizes the fifth dimension and x l , i = 1,2,3, the three noncompact dimen¬ 
sions. Cosmological isotropy excludes dtdx* or dydx 1 terms, and homogeneity ensures n 
and b are independent of x. The t, y part of the background metric may then be taken to 
be conformally flat. One can write the metric for two-dimensional Minkowski spacetime in 
Milne form so that the branes are located at y = ±y 0 and collide at t — 0 . By expressing 
the metric in locally Minkowski coordinates, T = t cosh y and Y = t sinhy, one sees that the 
collision speed is {V/c) = tanh 2y 0 and the relative rapidity of the collision is 2yo. As long 
as the bulk metric is regular at the brane collision and possesses cosmological symmetry, 
the line element may always be put into the form ©• Furthermore, by suitably rescaling 
coordinates one can choose 6(0,?/) = n(0,y) = 1. 

In order to describe perturbations about this background, one needs to specify an 
propriate gauge choice. Five-dimensional longitudinal gauge is particularly convenient 
firstly, it is completely gauge-fixed; secondly, the brane trajectories are unperturbed in this 
gauge { 32 )], so that the Israel matching conditions are relatively simple; and finally, in the 
absence of anisotropic stresses, the traceless part of the Einstein G* (spatial) equation yields 
a constraint amongst the perturbation variables, reducing them from four to three. In light 
of these advantages, we will work in five-dimensional longitudinal gauge for the entirety of 
this paper. 

Our three solution methods are as follows: 
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• Series expansion in t 

The simplest solution method for the background is to solve for the metric functions 
n(t, y ) and b(t, y ) as a series in powers of t about t — 0. At each order, the bulk Einstein 
equations yield a set of ordinary differential equations in y, with the boundary conditions 
provided by the Israel matching conditions. These are straightforwardly solved. A similar 
series approach, involving powers of t and powers of t times In t suffices for the perturbations. 

The series approach is useful at small times (t/L) <C 1 since it provides the precise solution 
for the background plus generic perturbations, close to the brane collision, for all y and for 
any collision rapidity ijq. It allows one to uniquely specify four-dimensional asymptotic data 
as t tends to zero. However, the series thus obtained fails to converge at quite modest times. 
Following the system to long times requires a more sophisticated method. Instead of taking 
(t/L) as our expansion parameter, we want to use the dimensionless rapidity of the brane 
collision y 0 , and solve at each order in y 0 . 

• Expansion in Dirichlet/Neumann polynomials in y 

In this approach we represent the spacetime metric in terms of variables obeying either 
Dirichlet or Neumann boundary conditions on the branes. We then express these variables 
as series of Dirichlet or Neumann polynomials in y and yo, bounded at each subsequent 
order by an increasing power of the collision rapidity y 0 . (Recall that the range of the y 
coordinate is bounded by \y\ < yo). The coefficients in these expansions are undetermined 
functions of t. By solving the five-dimensional Einstein equations perturbatively in yo, we 
obtain a series of ordinary differential equations in t, which can then be solved exactly. In 
this Dirichlet/Neumann polynomial expansion, the Israel boundary conditions on the branes 
are satisfied automatically at every order in yo, while the initial data at small t are provided 
by the previous series solution method. 

The Dirichlet/Neumann polynomial expansion method yields simple, explicit solutions 
for the background and perturbations as long as (t/L) is smaller than l/yo- Since yo <C 1, 
this considerably improves upon the naive series expansion in t. However, for (t/L) of order 
l/yo, the expansion fails because the growth in the coefficients overwhelms the extra powers 
of yo at successive orders. Since (t/L) ~ 1 /2/ 0 corresponds to brane separations of order 
the AdS radius, the Dirichlet/Neumann polynomial expansion method fails to describe the 
late-time behavior of the system, and a third method is needed. 


7 



• Expansion about the scaling solution 


The idea of our third method is to start by identifying a scaling solution, whose form is 
independent of y 0 for all y 0 • Cl. This scaling solution is well-behaved for all times and 
therefore a perturbation expansion in y 0 about this solution is similarly well-behaved, even 
at very late times. To find the scaling solution, we first change variables from t and y 
to an equivalent set of dimensionless variables. The characteristic velocity of the system 
is the brane speed at the collision, V = c tanh 2y 0 ~ 2q/o, for small y 0 , where we have 
temporarily restored the speed of light c. Thus we have the dimensionless time parameter 
x = y 0 ct/L Vt/L, of order the time for the branes to separate by one AdS radius. We 
also rescale the //-coordinate by defining u = y/yo, whose range is — 1 < uj < 1, independent 
of the characteristic velocity. 

As we shall show, when re-expressed in these variables, for small yo, the bulk Einstein 
equations become perturbatively ultralocal: at each order in y 0 one only has to solve an 
ordinary differential equation in u ;, with a source term determined by time derivatives of 
lower order terms. The original partial differential equations reduce to an infinite series of 
ordinary differential equations in u which are then easily solved order by order in yo- 

This method, an expansion in yo about the scaling solution, is the most powerful and may 
be extended to arbitrarily long times t and for all brane separations. As is well known for this 
model, a Birkhoff-like theorem holds for backgrounds with cosmological symmetry. The bulk 
in between the two branes is just a slice of five-dimensional Schwarzschild-AdS spacetime 


26 


27j within which the two branes move, with a virtual black hole lying outside of 


the physical region, beyond the negative-tension brane. As time proceeds, the negative- 
tension brane becomes closer and closer to the horizon of the Schwarzschild-AdS black hole. 
Even though its location in the Birkhoff-frame (static) coordinates freezes (see Figure 1), 
its proper speed grows and the y 0 expansion fails. Nonetheless, by analytic continuation 
of our solution in uj and x, we are able to circumvent this temporary breakdown of the yo 
expansion and follow the positive-tension brane, and the perturbations localized near it, as 
they run off to the boundary of anti-de Sitter spacetime. 

Our expansion about the scaling solution is closely related to derivative-expansion tech¬ 


niques developed earlier by a number of authors 


7,0,0- 


In these works, an expansion in 


terms of brane curvature over bulk curvature was used. For cosmological solutions, this is 





T 


1 b 

Figure 1: The background brane scale factors b± plotted as a function of the Birkhoff-frame 
time T, where b± have been normalized to unity at T = 0. In these coordinates the bulk is 
Schwarzschild-AdS: the brane trajectories are then determined by integrating the Israel matching 
conditions, given in Appendix E. In the limit as T —> oo, the negative-tension brane asymptotes to 
the event horizon of the black hole, while the positive-tension brane asymptotes to the boundary 
of AdS. 

equivalent to an expansion in L7i + , where 7i + is the Hubble constant on the positive-tension 
brane. However, we specifically want to study the time-dependence of the perturbations for 
all times from the narrowly-separated to the well-separated brane limit. For this purpose it is 
better to use a time-independent expansion parameter ( y 0 ) and to include all the appropriate 
time-dependence order by order in the expansion. 

Moreover, in the earlier works, the goal was to find the four-dimensional effective de¬ 
scription more generally, without specifying that the branes emerged from a collision with 
perturbations in the lowest Kaluza-Klcin modes. Consequently, the solutions obtained con¬ 
tained a number of undetermined functions. In the present context, however, the initial 
conditions along the extra dimension are completely specified close to the brane collision 
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by the requirement that only the lowest Kaluza-Klcin modes be excited. The solutions we 
obtain here are fully determined, with no arbitrary functions entering our results. 

Returning to the theme of the four-dimensional effective theory, we expect on general 
grounds that this should be valid in two particular limits: firstly, as we have already dis¬ 
cussed, a Kaluza-Klcin description will apply at early times near to the collision, when the 
separation of the branes is much less than L. Here, the warping of the bulk geometry and 
the brane tensions can be neglected. Secondly, when the branes are separated by many AdS 
lengths, one expects gravity to become localized on the positive-tension brane, which moves 
ever more slowly as time proceeds, so the four-dimensional effective theory should become 
more and more accurate. 

Equipped with our five-dimensional solution for the background and perturbations ob¬ 
tained by expanding about the scaling solution, we find ourselves able to test the four¬ 
dimensional effective theory explicitly. We will show that the four-dimensional effective the¬ 
ory accurately captures the five-dimensional dynamics to leading order in the ^-expansion, 
but fails at the Erst nontrivial order. Our calculations reveal that the four-dimensional 
perturbation modes undergo a mixing in the transition between the Kaluza-Klcin effective 
theory at early times and the brane-localized gravity at late times. This effect is a conse¬ 
quence of the momentary breakdown of the effective theory when the brane separation is 
of the order of an AdS length, and cannot be seen from four-dimensional effective theory 
calculations alone. 


III. SERIES SOLUTION IN t 


As described above, we End it simplest to work in coordinates in which the brane locations 
are fixed but the bulk evolves. The bulk metric is therefore given by (HD, with the brane 
locations fixed at y — ±yo for all time t. The five-dimensional solution then has to satisfy 
both the Einstein equations and the Israel matching conditions on the branes 0 - 

The bulk Einstein equations read G u = — AS”, where the bulk cosmological constant is 
A = —Q/L 2 (we work in units in which the four-dimensional gravitational coupling 87^4 = 
8ttG 5 /L = 1). Evaluating the linear combinations Gq + G\ and Gq + G 5 — (3/2 )G\ (where 0 
denotes time, 5 labels the y direction, and i denotes one of the noncompact directions, with 
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no sum implied), we find: 


v 


[3,tt ~ P,yy + /3 T — @ t y + 12 e 17 — 0 

1/ "° P\)-2e 2v = 0 , 


,TT 


z/ 


,w 4~ 2 (A 


( 2 ) 

(3) 


where (t/L) = eh (3 = 3 In 6 and z/ = In (nt/L). The Israel matching conditions on the 
branes read 31, -3J] 


(4) 


b ; y n tV nt 
b n L ’ 

where all quantities are to be evaluated at the brane locations, y = ±y 0 . 

We will begin our assault on the bulk geometry by constructing a series expansion in t 
about the collision, implementing the Israel matching conditions on the branes at each order 
in t. This series expansion in t is then exact in both y and the collision rapidity yo- It chief 
purpose will be to provide initial data for the more powerful solution methods that we will 
develop in the following sections. 

The Taylor series solution in t for the background was first presented in [3/ 


* 1 ■— -'*+<>(*) 


n = 1 + (sech y 0 sinh y) — + - sech 2/ 0 (—3 + 2 cosh 2 y + cosh 2 y 0 ) 

L 4 L z 


t i f 2 

6 = 1 + (sech y 0 sinh y) — -\— sech 2 y 0 (cosh 2 y — cosh 2 y 0 ) — + O 

L 2 L 


-i 

i?) 


(5) 

( 6 ) 


(Note that in the limit as t —> 0 we correctly recover compactihed Milne spacetime). 

Here, however, we will need the perturbations as well. Working in five-dimensional longi¬ 
tudinal gauge for the reasons given in the previous section, the perturbed bulk metric takes 
the form (see Appendix A) 

ds 2 = n 2 (-(1 + 2$ L ) dt 2 - 2W l dtdy + t 2 (1 - 2T L ) dy 2 ) + b 2 (1 - 2T L ) dx 2 , (7) 


with T^ = ~ d !l being imposed by the five-dimensional traceless G l - equation. The Israel 

matching conditions at y = ±yo then read 

nt ~y~ „ nt 


d>L, y = r L L , 
Performing a series expansion, we find 


®L, y - 


W L = 0. 


( 8 ) 


^ B B sech y 0 sinh y / B Bk 2 1 . 

<h L = -- +- j -- + (A - - - — + -Bk 2 In \kt\ 
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d> L = 
W L = 


-B cosh 2y (—1 + 6 y 0 coth 2y 0 ) sech 2 y 0 — -B sech 2 y 0 sinh 2 yj + 0(t ) (9) 

B sech y 0 sinh y 


B \ 1 

+ ( 2A — — + - Bk 2 In \kt\ + -B cosh 2y sech 2 y 0 ) + 0(f) (10) 


-B sech 2 y 0 (y cosh 2y — y 0 cosh 2y 0 sinh 2y) f + O(f 2 ), 


( 11 ) 
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where we have set L — 1 (to restore L, simply replace t t/L and k —> fcL). Except for a 
few specific instances, we will now adopt this convention throughout the rest of this paper. 
The two arbitrary constants A and B (which may themselves be arbitrary functions of k) 
have been chosen so that, on the positive-tension brane, to leading order in y 0 , goes as 

$ L = A - | + O(y 0 ) + 0(k 2 ) + 0(t). (12) 

IV. EXPANSION IN DIRICHLET/NEUMANN POLYNOMIALS 

A. Background 

Having solved the relevant five-dimensional Einstein equations as a series expansion in 
the time t before or after the collision event, we now have an accurate description of the 
behavior of the bulk at small t for arbitrary collision rapidities. However, in order to match 
onto the incoming and outgoing states, we really want to study the long-time behavior 
of the solutions, as the branes become widely separated. Ultimately this will enable us 
to successfully map the system onto an appropriate four-dimensional effective description. 
Instead of expanding in powers of the time, we approximate the five-dimensional solution 
as a power series in the rapidity of the collision, and determine each metric coefficient for 
all time at each order in the rapidity. 

Our main idea is to express the metric as a series of Dirichlet or Neumann polynomials 
in y 0 and y, bounded at order n by a constant times y% , such that the series satisfies the 
Israel matching conditions exactly at every order in y 0 . To implement this, we first change 
variables from b and n to those obeying Neumann boundary conditions. From Q. b/n is 
Neumann. Likewise, if we define N(t,y ) by 

1 

nt =-. 

N(t,y) -y 

then one can easily check that N(t,y ) is also Neumann on the branes. Notice that if N 
and b/n are constant, the metric (JTJ) is just that for anti-de Sitter spacetime. For fixed yo, 
N describes the the proper separation of the two branes, and b is an additional modulus 
describing the three-dimensional scale-factor of the branes. 

Since N and b/n obey Neumann boundary conditions on the branes, we can expand both 
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m a power series 


oo 


N = N 0 (t) + ^2 N n (t)Pn(y), b/n = q 0 (t)+ ^2q n (t)P n (y), (14) 


n =3 n =3 

where P n {y ) are polynomials 



satisfying Neumann boundary conditions and each bounded by \P n (y)\ < 2y^/(n — 2), for 


the relevant range of y. Note that the time-dependent coefficients in this ansatz may also be 
expanded as a power series in yo- By construction, our ansatz satisfies the Israel matching 
conditions exactly at each order in the expansion. The bulk Einstein equations are not 
satisfied exactly, but as the expansion is continued, the error terms are bounded by increasing 
powers of y 0 . 

Substituting the series ansatze (fill) into the background Einstein equations (j2j) and ([3fl , 
we may determine the solution order by order in the rapidity yo . At each order in yo, 
one generically obtains a number of linearly independent algebraic equations, and at most 
one ordinary differential equation in t. The solution of the latter introduces a number of 
arbitrary constants of integration into the solution. 

To fix the arbitrary constants, one first applies the remaining Einstein equations, allowing 
a small number to be eliminated. The rest are then determined using the series expansion 
in t presented in the previous section: as this solution is exact to all orders in y 0 , we need 
only to expand it out to the relevant order in yo, before comparing it term by term with 
our Dirichlet/Neumann polynomial expansion (which is exact in t but perturbative in y 0 ), 
taken to a corresponding order in t. The arbitrary constants are then chosen so as to ensure 
the equivalence of the two expansions in the region where both t and yo are small. This 
procedure suffices to fix all the remaining arbitrary constants. 

The Erst few terms of the solution are 



(16) 


(17) 


and 



93 = —2i 3 !/o + ..., 
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(18) 

(19) 



and the full solution up to 0{yl°) may be found in Appendix B. 


B. Perturbations 

Following the same principles used in our treatment of the background, we construct the 
two linear combinations 

04 = \{®l + *l), £4 = b 2 (y L - $ L ) = b 2 r L , (20) 

both of which obey Neumann boundary conditions on the branes, as may be checked from 
(j4j) and ( 0 ). In addition, Wl already obeys simple Dirichlet boundary conditions. 

The two Neumann variables, 04 and £ 4 , are then expanded in a series of Neumann poly¬ 
nomials and Wl is expanded in a series of Dirichlet polynomials, 

D n {y) = V n - Vo, n = 2,4,..., D n (y) = yD n _ 1 (y), n = 3,5,..., (21) 

each bounded by \D n (y)\ < y q for n even and r/g(n — 1 )/n n ^ n ~ 1) for n odd, over the relevant 
range of y. As in the case of the background, the time-dependent coefficients multiplying 
each of the polynomials should themselves be expanded in powers of 7 / 0 ■ 

To solve for the perturbations it is sufficient to use only three of the perturbed Einstein 
equations (any solution obtained may then be verified against the remainder). Setting 


< h L = 0 e “ 2l '- /3/3 ( 22 ) 

= -0 e - ^ 3 (23) 

W L = « r - 2 ^ 3 (24) 

where t = e r , 0 = 3In6 and v = hint, the Gf, G® and G\ equations take the form 

w )T = 2 0 iy — 4 e 3l//2 (0 e u/2 ) t y (25) 

0,r = w ,y ■* e3u ("0 e ~ V ),T (26) 

(tr^ /3 ),r = + (5 fP T - 5 /??„ - k 3 e^'A 

_ 2 e _ 2 „ +(J/3 ^^ T+p ^ ) _ wfjT ^y ( 27 ) 


Using our Neumann and Dirichlet ansatze for 04 , £4 and Wl, the Israel matching condi¬ 
tions are automatically satisfied and it remains only to solve (BSD, CU) and (EZD order by 
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order in the rapidity. The time-dependent coefficients for 0 4 , £4 and Wl are then found to 
obey simple ordinary differential equations, with solutions comprising Bessel functions in 
kt , given in Appendix C. Note that it is not necessary for the set of Neumann or Dirichlet 
polynomials we have used to be orthogonal to each other: linear independence is perfectly 
sufficient to determine all the time-dependent coefficients order by order in y 0 . 

As in the case of the background, the arbitrary constants of integration remaining in the 
solution after the application of the remaining Einstein equations are fixed by performing a 
series expansion of the solution in t. This expansion can be compared term by term with 
the series expansion in t given previously, after this latter series has itself been expanded 
in y 0 . The arbitrary constants are then chosen so that the two expansions coincide in the 
region where both t and // 0 are small. The results of these calculations, at long wavelengths, 
are: 


= A - B fk _ A-ln|fct j + ^At + j) y + 

Tl = 2A + B^\n\kt\ - (m + /j y - 
W L = ii-4 i 1 (jr 2 - j/J) + .. . 


(28) 

(29) 

(30) 


where the constants A and B can be arbitrary functions of k. The solutions for all k, to fifth 
order in y 0 , are given in Appendix C. 


V. EXPANSION ABOUT THE SCALING SOLUTION 

It is illuminating to recast the results of the preceding sections in terms of a set of 
dimensionless variables. Using the relative velocity of the branes at the moment of collision, 
V = 2c tanh y 0 — 2q/o (where we have temporarily re-introduced the speed of light c), we 
may construct the dimensionless time parameter x = yod/L ~ Vt/L and the dimensionless 
^/-coordinate u> — y/y 0 y(c/v). 

Starting from the full Dirichlet/Neumann polynomial expansion for the background given 
in Appendix B, restoring c to unity and setting t = xL/y 0 and y = ouy 0 , we find that 

n -1 = N(x) — ux + 0(yl) 

/ = q(x)+0(yl); 
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(31) 

(32) 



where 


N(x) 

q(x) 


1 - 
1 - 


x 2 3 x 4 25 x 6 343 x 8 
T ~ ~8 48 384 

3 x 2 7 x 4 55 x 6 245 x' 

~2 8 48 128" 


2187x 10 
1280 
4617 x 
~ 1280 


+ 0(x 12 ) 
o 

- + o(i 12 ). 


(33) 

(34) 


The single term in m linear in uj is necessary in order that n 1 satishes the correct boundary 
conditions. Apart from this one term, however, we see that to lowest order in yo the metric 
functions above turn out to be completely independent of uj. Similar results are additionally 
found for the perturbations. 

Later, we will see how this behavior leads to the emergence of a four-dimensional effective 
theory. For now, the key point to notice is that this series expansion converges only for 
x <C 1, corresponding to times t <C L/y 0 . In order to study the behavior of the theory 
for all times therefore, we require a means of effectively resumming the above perturbation 
expansion to all orders in x. Remarkably, we will be able to accomplish just this. The 
remainder of this section, divided into five parts, details our method and results: first, we 
explain how to find and expand about the scaling solution, considering only the background 
for simplicity. We then analyze various aspects of the background scaling solution, namely, 
the brane geometry and the analytic continuation required to go to late times, before moving 
on to discuss higher-order terms in the expansion. Finally, we extend our treatment to cover 
the perturbations. 


A. Scaling solution for the background 

The key to the our method is the observation that the approximation of small collision 
rapidity (y 0 <C 1) leads to a set of equations that are perturbatively ultralocal: transforming 
to the dimensionless coordinates x and u, the Einstein equations for the background © and 
a become 

+ P% ~ 12 e 2u = yl (x(x( 3 iX ) jX + x 2 /3 2 x ) (35) 

Yu; - ^ P% + 2e 2u = 7/o(x(xYY - ^ x 2 (5 2 x ) , (36) 

where we have introduced u — u + lny 0 . Strikingly, all the terms involving x-derivatives are 
now suppressed by a factor of y$ relative to the remaining terms. This segregation of x- and 
o;-derivatives has profound consequences: when solving perturbatively in yo, the Einstein 
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equations ( 1351 ) and © reduce to a series of ordinary differential equations in u, as opposed 
to the partial differential equations we started off with. 

To see this, consider expanding out both the Einstein equations ( 1351 ) and © as well 
as the metric functions (3 and i> as a series in positive powers of y 0 . At zeroth order in y 0 , 
the right-hand sides of © and ( 1351 ) vanish, and the left-hand sides can be integrated with 
respect to u to yield anti-de Sitter space. (This was our reason for using v = v + liiy 0 
rather than v\ the former serves to pull the necessary exponential term deriving from the 
cosmological constant down to zeroth order in y 0 , yielding anti-de Sitter space as a solution at 
leading order. As we are merely adding a constant, the derivatives of v and v are identical.) 
The Israel matching conditions on the branes ED, which in these coordinates read 




(37) 


are not however sufficient to fix all the arbitrary functions of x arising in the integration 
with respect to u. In fact, two arbitrary functions of x remain in the solution, which may 
be regarded as time-dependent moduli describing the three-dimensional scale factor of the 
branes and their proper separation. These moduli may be determined with the help of the 
Gg Einstein equation as we will demonstrate shortly. 

Returning to ( 1351 ) and ( 1351 ) at ?/q order now, the left-hand sides amount to ordinary 
differential equations in uj for the y^ corrections to (3 and v. The right-hand sides can no 
longer be neglected, but, because of the overall factor of y%, only the time-derivatives of (3 
and v at zeroth order in y 0 are involved. Since (3 and h have already been determined to 
this order, the right-hand sides therefore act merely as known source terms. Solving these 
ordinary differential equations then introduces two further arbitrary functions of x; these 
serve as y g corrections to the time-dependent moduli and may be fixed in the same manner 
as previously. 

Our integration scheme therefore proceeds at each order in y 0 via a two-step process: 
first, we integrate the Einstein equations (© and © to determine the ^-dependence of 
the bulk geometry, and then secondly, we fix the x-dependent moduli pertaining to the 
brane geometry using the Gf equation. This latter step works as follows: evaluating the Gf 
equation on the branes, we can use the Israel matching conditions © to replace the single 
^-derivatives that appear in this equation, yielding an ordinary differential equation in time 
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for the geometry on each brane. Explicitly, we find 

= 0, 


J ,x 


n 


( 38 ) 


,x 


b t± 1 

—— = ±—tanh y 0 , 
V o Vo 


(39) 


where five-dimensional considerations (see Section VI) further allow us to fix the constants 

of integration on the (±) brane as 

bb )X bb >t 
n y 0 n 

where the brane conformal time t± is defined on the branes via ndt = bdt±. When augmented 
with the initial conditions that n and b both tend to unity as x tends to zero (so that 
we recover compactified Milne spacetime near the collision), these two equations are fully 
sufficient to determine the two ^-dependent moduli to all orders in y 0 . 

Putting the above into practice, for convenience we will work with the Neumann variables 
N and q, generalizing (13T1) and (1321) to 

b 


n 1 = N(x,u) — lux, — = q(x,u). 

n 

Seeking an expansion of the form 

N(x,u) = No(x,uj) + yoNi(x,u)+ O(yo) 
q(x,u) = q 0 (x,u)+y%qi{x,u)+0(y%), 


(40) 

(41) 

(42) 


the Einstein equations (1331) and dD when expanded to zeroth order in yo immediately 
restrict N 0 and qo to be functions of x alone. The bulk geometry to this order is then simply 
anti-de Sitter space with time-varying moduli, consistent with m and (E2D- The moduli 
N 0 (x) and qo(x) may be found by integrating the brane equations (1331) . also expanded to 
lowest order in y 0 . In terms of the Lambert W-function 


W(x)e w< - x) = x, 


jd^ ]. W{x), defined implicitly by 

(43) 


the solution is 


N 0 (x) = e^ w{ ~ x2 \ q 0 (x ) = (1 + W(-x 2 )) e^~ x2 \ 


(44) 


Thus we have found the scaling solution for the background, whose form is independent of 
I/O, holding for any y 0 -C 1. Using the series expansion for the Lambert W-function about 
x = W{x) = 0, namely 3fj] 

(45) 


W(x) = J2 

m= 1 


-X 


m\ 
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W(x) 



Figure 2: The real values of the Lambert W-function. The solid line indicates the principal 
solution branch, Wq(x), while the dashed line depicts the W-\{x) branch. The two branches join 
smoothly at x = — 1/e where W attains its negative maximum of —1. 

we can immediately check that the expansion of our solution is in exact agreement with 
© and Ijm At leading order in y 0 then, we have succeeded in resumming the Dirich- 
let/Neumann polynomial expansion results for the background to all orders in x. 

Later, we will return to evaluate the corrections in our expansion about the scaling 
solution. In the next two subsections, however, we will first examine the scaling solution in 
greater detail. 

B. Evolution of the brane scale factors 

Using the scaling solution (©) to evaluate the scale factors on both branes, we find to 

O(y 2 0 ) 

b± = 1 ± xe~^ w{ ~ x2) = 1 ± \J—W(—x 2 ). (46) 

To follow the evolution of the brane scale factors, it is helpful to first understand the behavior 
of the Lambert W-function, the real values of which are displayed in Figure 2. For positive 
arguments the Lambert W-function is single-valued, however, for the negative arguments of 
interest here, we see that there are in fact two different real solution branches. The first 
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branch, denoted Wq(x), satisfies Wo(x) > — 1 and is usually referred to as the principal 
branch, while the second branch, W- 4 (x), is defined in the range W-\{x) < —1. The two 
solution branches join smoothly at x — —1/e, where W = —1. 

Starting at the brane collision where x = 0, the brane scale factors are chosen to satisfy 
b± = 1, and so we must begin on the principal branch of the Lambert W-function for 
which ILo(O) = 0. Thereafter, as illustrated in Figure 3, b + increases and decreases 
monotonically until at the critical time x = x c , when Wo(—x 2 ) = —1 and b_ shrinks to zero. 
From flUD, the critical time is therefore x c = e 2 = 0.606..., and corresponds physically to 
the time at which the negative-tension brane encounters the bulk black hole [3fi|. 

At this moment, the scale factor on the positive-tension brane has only attained a value 
of two. From the Birkhoff-frame solution, in which the bulk is Schwarzschild-AdS and the 
branes are moving, we know that the positive-tension brane is unaffected by the disappear¬ 
ance of the negative-tension brane and simply continues its journey out to the boundary of 
AdS. To reconcile this behavior with our solution in brane-static coordinates, it is helpful 
to pass to 1. 1 _, the conformal time on the positive-tension brane. Working to zeroth order in 
l/o, this may be converted into the dimensionless form 


x 4= Wt ±= y^ f ydt= I = xe-^~ x2) = sZ-Wi-x 2 ). (47) 

L L J b J q 0 (x) 

1 2 

Inverting this expression, we find that the bulk time parameter x = ^e - ^ 4 . The bulk 
time x is thus double-valued when expressed as a function of x 4 , the conformal time on the 
positive-tension brane: to continue forward in x 4 beyond x 4 = 1 (where x = x c ), the bulk 
time x must reverse direction and decrease towards zero. The metric functions, expressed 
in terms of x, must then continue back along the other branch of the Lambert W-function, 
namely the W-\ branch. I 11 this manner we see that the solution for the scale factor on 
the positive-tension brane, when continued on to the ILL 1 branch, tends to infinity as the 
bulk time x is reduced back towards zero (see dotted line in Figure 3), corresponding to the 
positive-tension brane approaching the boundary of AdS as x 4 —> 00 . 

For simplicity, in the remainder of this paper we will work directly with the brane confor¬ 
mal time x 4 itself. With this choice, the brane scale factors to zeroth order in yo are simply 
b± — 1 ± x 4 . 
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Figure 3: The scale factors b± on the positive-tension brane (rising curve) and negative-tension 
brane (falling curve) as a function of the bulk time parameter x, to zeroth order in r/o- The 
continuation of the positive-tension brane scale factor on to the W-\ branch of the Lambert W- 
function is indicated by the dashed line. 

C. Analytic continuation of the bulk geometry 


In terms of x±, the metric functions n and b are given by 


n = 


I t 2 

g 2 4 
1 — XX 4 


+ 0 {y j 


ob 


6 = ^ 


1 — XX ^ 


+ 0(y, 


(48) 


At X 4 = 1, the three-dimensional scale factor b shrinks to zero at all values of x except 
x — 1 (he. the positive-tension brane). Since b is a coordinate scalar under transformations 
of X 4 and x, one might be concerned that that the scaling solution becomes singular at 
this point. However, when we compute the y ( 2 corrections as we will do shortly, we will 
find that the corrections become large close to X\ = 1, precipitating a breakdown of the 
small- 7/0 expansion. Since it will later turn out that the scaling solution maps directly on to 
the four-dimensional effective theory, and that this, like the metric on the positive-tension 
brane, is completely regular at X 4 = 1, we are encouraged to simply analytically continue 
the scaling solution to times £4 > 1 . 

When implementing this analytic continuation careful attention must be paid to the 
range of the coordinate ix. Thus far, for times £4 < 1, we have regarded w as a coordinate 
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X 4 

(0 = 3 w = 2 o)=1.5 a) 



= 1.1 

= 1 


b 


Figure 4: The contours of constant uj in the ( b , X 4 ) plane. Working to zeroth order in yo, these 
are given by X4 = ^ ( buj ± \Jb 2 uj 2 — 4(6 — 1 )), where we have plotted the positive root using a solid 
line and the negative root using a dashed line. The negative-tension brane is located at uj = — 1 
for times X 4 < 1 , and the trajectory of the positive-tension brane is given (for all time) by the 
positive root solution for u = 1. The region delimited by the trajectories of the branes (shaded) 
then corresponds to the bulk. From the plot we see that, for 0 < X 4 < 1, the bulk is parameterized 
by values of uj in the range — 1 < uj < 1. In contrast, for X 4 > 1, the bulk is parameterized by 
values of uj in the range u> > 1 . 

spanning the fifth dimension, taking values in the range The two metric 

functions n and b were then expressed in terms of the coordinates X4 and uj. Strictly 
speaking, however, this parameterization is redundant: we could have chosen to eliminate u 
by promoting the three-dimensional scale factor b from a metric function to an independent 
coordinate parameterizing the fifth dimension. Thus we would have only one metric function 
n, expressed in terms of the coordinates X4 and b. 

While this latter parameterization is more succinct, its disadvantage is that the locations 
of the branes are no longer explicit, since the value of the scale factor b on the branes is time- 
dependent. In fact, to track the location of the branes we must re-introduce the function 
uj(x4,b) = {b + x 4 — l)/bx4 (inverting (|48l) at lowest order in y 0 ). The trajectories of the 
branes are then given by the contours u = ±1. 
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Figure 5: The three-dimensional scale factor 6 , plotted to zeroth order in yo as a function of £4 
and u), for £4 < 1 (left) and £4 > 1 (right). The positive-tension brane is fixed at w = 1 for 
all time (note the evolution of its scale factor is smooth and continuous), and for £4 < 1, the 
negative-tension brane is located at u = — 1 . 

The contours of constant u as a function of £4 and b are plotted in Figure 4. The analytic 
continuation to times £4 > 1 has been implemented, and the extent of the bulk is indicated 
by the shaded region. From the figure, we see that, if we were to revert to our original 
parameterization of the bulk in terms of £4 and u, the range of u required depends on 
the time coordinate £ 4 : for early times X 4 < 1 , we require only values of u in the range 
— 1 < uj < 1, whereas for late times £4 > 1, we require values in the range tv > 1. Thus, while 
the positive-tension brane remains fixed at u = 1 throughout, at early times £4 < 1 the value 
of u decreases as we head away from the positive-tension brane along the fifth dimension, 
whereas at late times £4 > 1, the value of u increases away from the positive-tension brane. 

While this behavior initially appears paradoxical if u is regarded as a coordinate along 
the fifth dimension, we stress that the only variables with meaningful physical content 
are the brane conformal time £4 and the three-dimensional scale factor b. These physical 
variables behave sensibly under analytic continuation. In contrast, u is simply a convenient 
parameterization introduced to follow the brane trajectories, with the awkward feature that 
its range alters under the analytic continuation at £4 = 1 . 

For the rest of this paper, we will find it easiest to continue parameterizing the bulk in 
terms of £4 and u, adjusting the range of the u where required. Figure 5 illustrates this 
approach: at early times £4 < 1 the three-dimensional scale factor b is plotted for values of 
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u in the range —1 <u<1 . At late times x\ > 1, we must however plot b for values of u in 
the range u > 1. In this fashion, the three-dimensional scale factor b always decreases along 
the fifth dimension away from the brane. 

We have argued that the scaling solution for the background, obtained at lowest order in 
y 0 , may be analytically continued across x^ — 1. There is a coordinate singularity in the x 4 , u 
coordinates but this does not affect the metric on the positive-tension brane which remains 
regular throughout. The same features will be true when we solve for the cosmological 
perturbations. The fact that the continuation is regular on the positive-tension brane and 
precisely agrees with the predictions of the four-dimensional effective theory provides strong 
evidence for its correctness. Once the form of the the background and the perturbations 
have been determined to lowest order in y 0 , the higher-order corrections are obtained from 
differential equations in y with source terms depending only on the lowest order solutions. 
It is straightforward to obtain these corrections for x± < 1. If we analytically continue them 
as described to x± > 1 as described, we automatically solve the bulk Einstein equations 
and the Israel matching conditions on the positive tension brane for all x 4 . The continued 
solution is well-behaved after the collision in the vicinity of the positive-tension brane, out 
to large distances where the y 0 expansion fails. 

D. Higher-order corrections 

In this section we explicitly compute the y% corrections. The size of these corrections 
indicates the validity of the expansion about the scaling solution, which perforce is only 
valid when the y\ corrections are small. 

Following the procedure outlined previously, we first evaluate the Einstein equations 
m> and m to 0 (y q) using the ansatze (BT1) and (TE21) . along with the solutions for No(x) 
and qo(x) given in (jUJ). The result is two second-order ordinary differential equations in u, 
which may straightforwardly be integrated yielding N { (x, uj) and q±(x, u) up to two arbitrary 
functions of x 4 . These time-dependent moduli are then fixed using the brane equations 4HD, 
evaluated at 0 (?/q) higher than previously. 
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Figure 6: The ratio of the y q corrections to the leading term in the small-yo expansion for b, 
plotted for X 4 < 1 (left) and X 4 > 1 (right), for the case where yo = 0.1. Where this ratio becomes 
of order unity the expansion about the scaling solution breaks down. The analogous plots for n 
display similar behavior. 


We obtain the result: 

W 2 It 2 2 

e2 *4 £2 4 yo 

17,1x4,0;) = -1-o- j 

I — UJX 4 30(—1 + a;x 4 )"(—1 + x|) 

( 0 : 4 ( 50 ; (-3 + a; 2 ) - 5 x 4 + 40a; (-3 + a; 2 ) x\ - 5 (-14 + 9a; 2 (-2 + a; 2 )) x^ 
+3o; 3 (-5 + 3ur) x\ - 19x 3 + 5 X 4 ) - 5(— 1 + x 2 ) 3 ln(l - x 2 )) + 0{y$), (49) 

, \ 1 -^4 ^'4?7o 

0(X4, to) — -1- 7 , - 

1 — 02 x 4 30(—l+o;x4) (—I+X 4 ) 

( - 5a; (-3 + a; 2 ) - 20x 4 + 5a; (-7 + 4a; 2 ) x\ - 10 (l - 12a ; 2 + 3a; 4 ) x^ 

+3a; (—20 — 5a ; 2 + 2a; 4 ) x 4 — 12x\ + 31a;x4 — 5a jx\ 

—5(—1 + x 2 ) 2 (a; - 2 x 4 + ujxfj ln(l - x 2 )) + 0{y q). (50) 


In Figure 6 , we have plotted the ratio of the y\ corrections to the corresponding terms 
at leading order: where this ratio becomes of order unity the expansion about the scaling 
solution breaks down. Inspection shows there are two such regions: the first is for times close 
to X 4 = 1 , for all ou, and the second occurs at late times X 4 > 1 , far away from the positive- 
tension brane. In neither case does the failure of the yo expansion indicate a singularity of 
the background metric: from the bulk-static coordinate system we know the exact solution 
for the background metric is simply Schwarzschild-AdS, which is regular everywhere. The 
exact bulk-static solution in Birkhoff frame tells us that the proper speed of the negative- 
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tension brane, relative to the static bulk, approaches the speed of light as it reaches the 
event horizon of the bulk black hole. It therefore seems plausible that a small-7/o expansion 
based upon slowly-moving branes must break down at this moment, when x 4 = 1 in our 
chosen coordinate system. 

Analytically continuing our solution in both x 4 and u around x 4 = 1 , the logarithmic 
terms in the t/q corrections now acquire imaginary pieces for > 1. However, since these 
imaginary terms are all suppressed by a factor of t/q, they can only enter the Einstein- 
brane equations (expanded as a series to t/q order) in a linear fashion. Hence the real and 
imaginary parts of the metric necessarily constitute independent solutions, permitting us to 
simply throw away the imaginary part and work with the real part alone. As a confirmation 
of this, it can be checked explicitly that replacing the In (1 — x±) terms in (1491) and (EDI) with 
In |1 — £ 4 1 still provides a valid solution to 0(t/q) of the complete Einstein-brane equations 
and boundary conditions. 

Finally, at late times where X 4 > 1, note that the extent to which we know the bulk 
geometry away from the positive-tension brane is limited by the t/q corrections, which become 
large at an increasingly large value of u, away from the positive-tension brane (see Figure 
6 ). The expansion about the scaling solution thus breaks down before we reach the horizon 
of the bulk black hole, which is located at u —> 00 for x 4 > 1. 

E. Treatment of the perturbations 

Having determined the background geometry to 0(t/q) in the preceding subsections, we 
now turn our attention to the perturbations. In this subsection we show how to evaluate the 
perturbations to 0(y q) by expanding about the scaling solution. The results will enable us 
to perform stringent checks of the four-dimensional effective theory and moreover to evaluate 
the mode-mixing between early and late times. 

In addition to the dimensionless variables x = y^ct/L and u = y/yo, when we consider the 
metric perturbations we must further introduce the dimensionless perturbation amplitude 
B = By 2 c 2 /L 2 ~ BV 2 /L 2 and the dimensionless wavevector k = kL/y 0 ~ ckL/V. In this 
fashion, to lowest order in 7/0 and k, we then find <f>£ = A — B/t 2 = A — B/x 2 and similarly 
kct = kx. (Note that the perturbation amplitude A is already dimensionless however). 

Following the treatment of the perturbations in the Dirichlet/Neumann polynomial ex- 
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pansion, we will again express the metric perturbations in terms of Wl (obeying Dirichlet 
boundary conditions), and the Neumann variables 04 and £4, defined in ( 1201 ) . Hence we seek 
an expansion of the form 


(p 4 (x 4 } uj) 

£ 4 {x 4 ,u) 

W L (x 4 ,u) 


040 {x 4 ,iv) +yl(j) 4 i(x 4 ,u ) + 0 ( 2 / 0 ), 
£40(2+ u) + y^ 41 (x 4 ,uo) + O(yo), 
W L 0 (x 4 ,u) + yl W L1 {x 4 ,u) + 0 ( 2 / 0 ). 


(51) 

(52) 

(53) 


As in the case of the background, we will use the G\ equation evaluated on the brane to fix 
the arbitrary functions of x 4 arising from integration of the Einstein equations with respect 
to to. By substituting the Israel matching conditions into the G | equation, along with the 
boundary conditions for the perturbations, it is possible to remove the single (^-derivatives 
that appear. We arrive at the following second-order ordinary differential equation, valid on 
both branes, 


0 = 2n{x\ — 1)(26 2 0 4 + £4 )b 2 + b 2 (nx 4 (x 2 — 3) — (a ; 2 — 1 )n) (2 & 2 £> 4 — £ 4 ) 

+bb ^4 nx 4 (xl - 3)(6 2 0 4 + £ 4 ) - (x 2 4 - 1)(4(6 2 0 4 + £ 4 )h - n(lO5 2 0 4 + £ 4 )) j 
+bn(xl - 1) ^4(5 2 0 4 + £4 )b + 25 3 0 4 - 5£ 4 j , (54) 


where dots indicate differentiation with respect to x 4 , and where, in the interests of brevity, 
we have omitted terms of 0 (k 2 ). 

Beginning our computation, the G\ and G\ Einstein equations when evaluated to lowest 
order in y 0 immediately restrict 040 and £40 to be functions of x 4 only. Integrating the 
G® equation with respect to to then gives Wlq in terms of 0 4O and £ 40 , up to an arbitrary 
function of x 4 . Requiring that H + 0 vanishes on both branes allows us to both fix this 
arbitrary function, and also to solve for £40 in terms of £>40 alone. Finally, evaluating (ED on 
both branes to lowest order in yo and solving simultaneously yields a second-order ordinary 
differential equation for £> 40 , with solution 


£>40 




B 

2 x 2 4 


+ 0 (k 2 ), 


(55) 


where the two arbitrary constants have been chosen to match the small-f series expansion 
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Figure 7: The perturbation to the three-dimensional scale factor, b 2 ^L, plotted on long wave¬ 
lengths to zeroth order in yo f° r early times (left) and late times (right). Only the B mode is 

displayed ( i.e. A = 0 and B = 1). Note how the perturbations are localized on the positive-tension 
brane (located at u = 1 ), and decay away from the brane. 

given in Section III. With this choice, 

„ 115 B ( „ 3 B\ 9 ~ 9n , . 

£40 — —A H—-- ^2 + ( ^ g - J x 4 + 0{k ) (56) 

W L0 = ^ (3Axl(-2 + oox 4 ) + Bx 4 (^X 4 + u(1 - + 0(k 2 ). (57) 

(1 - xlY \ 4 8 / 

The resulting behavior for the perturbation to the three-dimensional scale factor, b 2 x $> l, is 
plotted in Figure 7. 

In terms of the original Newtonian gauge variables, an identical calculation (working now 
to all orders in k) yields, 

= ^T^TT- i A o J o(kx 4 ) + BoYoCkx,)) 
o [x^ i) 

+ — + —:- 2 ^ 1 { A oJi(kx4) + (^4)) + 0(?/o), (58) 

X 4 y i x 4 / 

= 3 l _ ^2fc(l - ujx 4 ) 2 (A 0 J 0 (kxY + B 0 Y 0 (kx 4)) 

-3 (X 4 +u(—2 + ux 4 ))(A 0 Ji(kx 4 ) + B 0 Yi(kx 4 ))^j +0{yl), (59) 

W L = 2 x\ e~^ (fc (1 - ujx 4) (A 0 J 0 (kx 4 ) + B 0 Y 0 (kx 4)) 

+u; (A 0 Ji(kx 4 ) + 5 0 1i(^4))^ + 0(yjj), 
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( 60 ) 











where the constants Aq and Bq are given by 

A) = -f-—7- + -Bk(\n 2 — 7) + 0(2/0)) Bo = —- —fO(j/g). ( 61 ) 

k 8k * 4 

To evaluate the t/q corrections, we repeat the same sequence of steps: integrating the 
Gf and G\ Einstein equations (at yl higher order) gives us <f> 4 i and £41 up to two arbitrary 
functions of x 4 , and integrating the G® equation then gives us Wj A in terms of these two 
arbitrary functions plus one more. Two of the three arbitrary functions are then determined 
by imposing the Dirichlet boundary conditions on Wj A , and the third is found to satisfy a 
second-order ordinary differential equation after making use of AMD on both branes. This 
can be solved, and the constants of integration appearing in the solution are again chosen 
so as to match the small-f series expansion of Section III. 

Converting back to the original longitudinal gauge variables, the results to 0(y q) and to 
0(k 2 ) take the schematic form 

= fo + vl(ft + ft ln (! + x d + ft In (1 - x A ) + ft In (1 - ux 4 ), (62) 

= ff + yliff + ft In (1 + Z 4 ) + ft In (1 - x 4 ) + ft 111 (1 - ux 4 ), (63) 

W L = e“H (/ 0 M/ + y *(fW + ff ln (1 + x 4 ) + ff In (1 - x 4 ) + /]' ln (1 - ux 4 ))) , (64) 

where the / are rational functions of x 4 and u which, due to their length, have been listed 

separated separately in Appendix D. (If desired, more detailed results including the 0(k 2 ) 
corrections are available j30|). 

It is easy to check that the results obtained by expanding about the scaling solution are 
consistent with those obtained using our previous method based upon Dirichlet/Neumann 
polynomials. Taking the results from the polynomial expansion given in Appendix C, sub¬ 
stituting t = (x 4 /yo)e-2 x 4 and y = uyo, retaining only terms of 0{;yf) or less, one finds 
agreement with the results listed in Appendix D after these have been re-expressed as a se¬ 
ries in x 4 . This has been checked explicitly both for the background and the perturbations. 

Just as in the case of the background, the small-yo expansion breaks down for times close 
to x 4 — 1 when the yl corrections to the perturbations become larger than the corresponding 
zeroth order terms. Again, we will simply analytically continue the solution in x 4 and u 
around this point. In support of this, the induced metric on the positive-tension brane is, 
to zeroth order in y 0 , completely regular across x 4 — 1 , even including the perturbations as 
can be seen from © and ©. 
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As in the case of the background, any imaginary pieces acquired from analytically con¬ 
tinuing logarithmic terms are all suppressed by order y^. Thus they may only enter the 
Einstein-brane equations (when these are expanded to order y$) in a linear fashion, and 
hence the real and imaginary parts of the metric constitute independent solutions. We 
can therefore simply drop the imaginary parts, or equivalently replace the In (1 — x 4 ) and 
In (1 — tux 4 ) terms with In 11 — x 4 \ and In 11 — ujx 4 ) respectively. We have checked explicitly 
that this still satisfies the Einstein-brane equations and boundary conditions. 

VI. COMPARISON WITH THE FOUR-DIMENSIONAL EFFECTIVE THEORY 


We have now arrived at a vantage point from which we may scrutinize the predictions of 
the four-dimensional effective theory in light of our expansion of the bulk geometry about 
the scaling solution. We will find that the four-dimensional effective theory is in exact 
agreement with the scaling solution. Beyond this, the y$ corrections lead to effects that 
cannot be described within a four-dimensional effective framework. Nonetheless, the higher- 
order corrections are automatically small at very early and very late times, restoring the 
accuracy of the four-dimensional effective theory in these limits. 

In the near-static limit, the mapping from four to five dimensions may be calculated from 
the moduli space approach 0 , ItI [fisl ]: putting the four-dimensional effective theory metric 


into Einstein frame, the mapping reads 


9% = cosh 2 (<I>/V6)g%, g uv = sinh 2 ((j)/V6)g*„, 


1tu 




1 l±V> 


(65) 


where g+ u and g~ v are the metrics on the positive- and negative-tension branes respectively, 
and 0 is the radion. Two of us have shown elsewhere that on symmetry grounds this is the 
unique local mapping involving no derivatives and that to leading order the action for 
g^ u and 0 is that for Einstein gravity with a minimally coupled scalar held. 

Solving the four-dimensional effective theory is trivial: the background is conformally 
hat, g* v = 64 ( 14 ) and the Einstein-scalar equations yield the following solution, unique 
up to a sign choice for 00 , of the form 

64 = C 4 t 4 , eV^ 10 = A 4 t 4 , (66) 

with 0o the background scalar held, and A 4 and C 4 arbitrary constants. (Throughout this 
paper we adopt units where 87 tG 4 = 1 ). 
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According to the map m, the brane scale factors are then predicted to be 

1 . _^o 


b+ = 


-b 4 e |l i e\Z7^ 0 ^ = 1 ± ^ 4 ^ 4 , 


(67) 


where we have chosen Ca = 4A 4 , so that the brane scale factors are unity at the brane colli¬ 
sion. As emphasized in f22j, the result m is actually exact for the induced brane metrics, 
when t 4 is identihed with the conformal time on the branes. From this correspondence, one 
can read off the five-dimensional meaning of the parameter A 4 : it equals L~ l tanh y 0 (our 


32] by a factor of 2). 


definition of y 0 differs from that of t _ a 

With regard to the perturbations, in longitudinal gauge (see e.g. 
line element of the four-dimensional effective theory reads 


41]) the perturbed 


dsj = b'j(t 4 ) [-(1 + 2$ 4 )dt 4 + (1 - 2<& 4 )dx 2 ] , 


( 68 ) 


and the general solution to the perturbation equations at small k 


42 is 


$4 — — ^Ao J\(kt 4 ) + BoYi(kt 4 )j , (69) 

= 7^k (AoMkt*) + B 0 Y 0 (ktS) ~ J- (A 0 Ji(kt A ) + B 0 Y 1 (kt 4 )^ , (70) 

with Aq and Bq being the amplitudes of the two linearly independent perturbation modes. 


A. Background 

In the case of the background, we require only the result that the scale factors on the 
positive- and negative-tension branes are given by 

b± = 1 ± A 4 t 4 , (71) 

where the constant A 4 = L~ l tanhy 0 and t 4 denotes conformal time in the four-dimensional 
effective theory. (Note this solution has been normalized so as to set the brane scale factors at 
the collision to unity). Consequently, the four-dimensional effective theory restricts 6 + +6_ = 
2. In comparison, our results from the expansion about the scaling solution PI) give 

6+ + 6_ ^ 2 + -f) 3 y ° + 0(»i). (72) 

Thus the four-dimensional effective theory captures the behavior of the full theory only in 
the limit in which the corrections are small, i.e., when the scaling solution is an accurate 
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description of the higher-dimensional dynamics. At small times such that 14 C 1 , the y$ 
corrections will additionally be suppressed by 0(xl), and so the effective theory becomes 
increasingly accurate in the Kaluza-Klcin limit near to the collision. Close to x 4 = 1, the 
small-r/o expansion fails hence our results for the bulk geometry is no longer reliable. For 
times £4 > 1 , the negative-tension brane no longer exists and the above expression is not 
defined. 

We can also ask what the physical counterpart of £ 4 , conformal time in the four- 
dimensional effective theory, is: from CT, we find 

t A = — 2 -r b ~ = ^ - ——- 3 - 2 3 ( x 4(5 - 14 x\ + 5:4) - 5x 4 (—1 + x\f ln(l - x 4 )) + 0(yl). 

4 t/ 0 4 vJ ^ 4 PC /j^ j 

(73) 

In comparison, the physical conformal times on the positive- and negative-tension branes, 
defined via bdt± — ndt — (n/yo)(l — x'l) e~ x dx 4 , are, to 0(y q), 

t + = —4-—-tt flO — 30 x 4 — x\(h — ldx\ + 5xi) + 5xJlxj) 2 ln(l — x 4 )') (74) 

Vo 30(1 — xj) 3 ' > 

t_ = —-—- n (10 — 30:4 + xl(5 — 14:4 + 5 ^ 4 ) — 5x 4 (1 ■«- x 2 X ln(l — 4)1, (75) 

Vo 30(1 — xl) 3 ' ) 

where we have used (SI and m- 

Remarkably, to lowest order in r/o, the two brane conformal times are in agreement not 
only with each other, but also with the four-dimensional effective theory conformal time. 
Hence, in the limit in which y\ corrections are negligible, there exists a universal four¬ 
dimensional time. In this limit, t 4 = x^/y^ and the brane scale factors are simply given by 
b± = 1 ± A 4 f 4 = 1 ± x 4 . The four-dimensional effective scale factor, & 4 , is given by 

(& 4 ) 2 = b 2 + — b 2 _ = 4 A 4 f 4 = Ax\. (76) 

In order to describe the full five-dimensional geometry, one must specify the distance 
between the branes d as well as the metrics induced upon them. The distance between the 
branes is of particular interest in the cyclic scenario, where an inter-brane force depending 
on the inter-brane distance d is postulated. In the lowest approximation, where the branes 
are static, the four-dimensional effective theory predicts that 

d =i,ncoth (H) =ihi (r ) =Lin (444) - (77) 
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Substituting our scaling solution and evaluating to leading order in y 0l we find 




(78) 


(Again, this quantity is ill-defined for X 4 > 1). 

In the full five-dimensional setup, a number of different measures of the inter-brane dis¬ 
tance are conceivable, and the inter-brane force could depend upon each of these, according 
to the precise higher-dimensional physics. One option would be to take the metric distance 
along the extra dimension 

r y o rvo /■> , 2 

d m = LI y/g^dy = L / ntdy = L / nx^e~^ XA duj. (79) 

J—yo J ~yo J—l 

Using (HID, we obtain 


dm = L f 1 — —do; + 0(y 2 ) =Lln(\±^) + O(y 2 0 ), (80) 

7 -i J- (XX4 \ i X4 j 

in agreement with (ZED- 

An alternative measure of the inter-brane distance is provided by considering affinely 
parameterized spacclikc geodesics running from one brane to the other at constant Birkhoff- 
frame time and noncompact coordinates x l . The background interbrane distance is just the 
affine parameter distance along the geodesic, and the fluctuation in distance is obtained 
by integrating the metric fluctuations along the geodesic, as discussed in Appendix E. One 
finds that, to leading order in y 0 only, the geodesic trajectories lie purely in the ^-direction. 
Hence the affine distance d a is trivially equal to the metric distance d m . at leading order, 
since 

d a = L J a/ g a bX a x b dA = L J ntydX = d m , (81) 

where the dots denote differentiation with respect to the affine parameter A. Both measures 
of the inter-brane distance therefore coincide and are moreover in agreement with the four¬ 
dimensional effective theory prediction, but only at leading order in y 0 . 


B. Perturbations 

Since the four-dimensional Newtonian potential <f> 4 represents the anticonformal part of 
the perturbed four-dimensional effective metric (see (I68[) i. it is unaffected by the conformal 
factors in (ESI) relating the four-dimensional effective metric to the induced brane metrics. 
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Hence we can directly compare the anticonformal part of the perturbations of the induced 
metric on the branes, as calculated in five dimensions, with 2 <f >4 in the four-dimensional 
effective theory. The induced metric on the branes is given by 


ds 2 = b 2 (-(1 + 2*J.) di| + (1 - 2<S l ) dx 2 ) 

= b 2 (1 + <S> L - (-(1 + dt| + (1 - (f L + S t )) d2 2 ) , 


(82) 

(83) 


where the background brane conformal time, t±, is related to the bulk time via bdt± = ndt. 
The anticonformal part of the metric perturbation is thus simply <£>L + 'I r L- If is this quantity, 
evaluated on the branes to leading order in y 0 , that we expect to correspond to 2<3> 4 in the 
four-dimensional effective theory. 

Using our results m and m from expanding about the scaling solution, we have to 

O(y 2 0 ), 

+ T l ) + = + T l )_ = (j> 40 (x 4 ) = (H 0 Ji(fcr 4 ) + BoY^kx^ , (84) 

with Ho and Bq as given in (ED- On the other hand, the Newtonian potential of the 
four-dimensional effective theory is given by (EUD- Since £ 4 , the conformal time in the four¬ 
dimensional effective theory, is related to the physical dimensionless brane conformal time 
x 4 by £ 4 = x 4 /y Q (to lowest order in y 0 ), and moreover k = k/y 0 , we have kx 4 = kt 4 . 
Hence the four-dimensional effective theory prediction for the Newtonian potential is in 
exact agreement with the scaling solution holding at leading order in y 0l upon identifying 
Ho with Ho/ 2/0 and Bq with B^/y^. The behavior of the Newtonian potential is illustrated 
in Figure 8. 

Turning our attention now to the radion perturbation, <50, we know from our earlier con¬ 
siderations that this quantity is related to the perturbation Sd in the inter-brane separation. 
Specifically, from varying (ED, we find 


5d = 2 L cosech 

Inserting the four-dimensional effective theory predictions for 0 and <50, we obtain 

— = ^^ t )2 -^ (AoJo(kt 4 ) + B 0 Y 0 (kt 4 )^j — — ^ A 0 Ji(kt 4 ) + B 0 Yi(kt 4 ))^j , ( 86 ) 

where to lowest order H 4 = y Q + 0(yl). 
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Figure 8 : The four-dimensional Newtonian potential $4 on the positive-tension brane, plotted to 
zeroth order in 7/0 as a function of the time x 4 for wavelength k = 1. The plot on the left illustrates 
the mode with A = 1 and B = 0, while the plot on the right has A = 0 and B = 1. 


In comparison, the perturbation in the metric distance between the branes is 


5 d r , 


cyo 


ntr L dy = 


' -y 0 


'-1 


1 nx 4 U _i x 2 

— e 1 ‘dw, 


(87) 


where we have used 02DD- Evaluating the integral using m, to an accuracy of 0(y g) we 
obtain 


5d m _ f 1 nx 4 $, 40 (x 4 ) _i I 2 J _ 2x 4 ^4o(x 4 ) 

L ~ y_4 6 2 e 2 (1-x 2 ) 2 

= ^2Y) ^ kx4,(A 0 J 0 (kx 4 ) + B 0 Y 0 (kx 4 )) - 4(A 0 Ji(A:x4 ) + B 0 Y\(kx 4 ))^ , (88) 

which is in agreement with m when we set ^ 4^4 rv./ yoU = x 4 , along with A 0 = A 0 /y 0 , B 0 = 
B 0 /y 0 and k = ky 0 . The calculations in Appendix E show moreover that the perturbation 
in the affine distance between the branes, 5d a , is identical to the perturbation in the metric 
distance 5d m , to lowest order in y 0 . 

The four-dimensional effective theory thus correctly predicts the Newtonian potential 
$4 and the radion perturbation 5<fr, but only in the limit in which the corrections are 
negligible and the bulk geometry is described by the scaling solution. While these corrections 
are automatically small at very early or very late times, at intermediate times they cannot be 
ignored and introduce effects that cannot be described by four-dimensional effective theory. 
The only five-dimensional longitudinal gauge metric perturbation we have not used in any 
of the above is Wi : this component is effectively invisible to the four-dimensional effective 
theory, since it vanishes on both branes and has no effect on the inter-brane separation. 
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VII. MIXING OF GROWING AND DECAYING MODES 


Regardless of the rapidity of the brane collision y 0 , one expects a four-dimensional effective 
description to hold both near to the collision, when the brane separation is much less than 
an AdS length, and also when the branes are widely separated over many AdS lengths. In 
the former case, the warping of the bulk geometry is negligible and a Kaluza-Klein type 
reduction is feasible, and in the latter case, one expects to recover brane-localized gravity. 
At the transition between these two regions, when the brane separation is of order one AdS 
length, one might anticipate a breakdown of the four-dimensional effective description. 

However, when the brane separation is of order a few AdS lengths, the negative-tension 
brane reaches the horizon of the bulk black hole and the small-y 0 expansion fails. This 
failure hampers any efforts to probe the breakdown of the four-dimensional effective theory 
at x A = 1 directly; instead, we will look for evidence of mixing between the four-dimensional 
perturbation modes in the transition from Kaluza-Klein to brane-localized gravity. 

To see this in action we have to compare the behavior of the perturbations at very 
small times with that at very late times: in both of these limits a four-dimensional effective 
description should apply, regardless of the collision rapidity y 0 , in which the four-dimensional 
Newtonian potential $4 satisfies 


$4 


£4 


A Q Ji{kx A ) + B 0 Y\(kx A ) 


(89) 


Expanding this out on long wavelengths k ■ Cl, taking in addition kx 4 < 1, we find 

$4 = - ^ B A k 2 In kx A - ^ A° 4 k 2 x\ + 0(fc 4 ), (90) 

where the dimensionless constants A 4 and B® are given in terms of the five-dimensional 
perturbation amplitudes A and B by 


A 0 = -A-—B- -Bk 2 , B° 4 = -B, 

4 2 16 8 4 2 ’ 


(91) 


where we have used m and (JED, recalling that $4 = 0 4O at leading order in y 0 . 

In comparison, using our results from expanding about the scaling solution, we find the 
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Newtonian potential on the positive-tension brane at small times X 4 -C 1 is given by 


$4 = - 


B + - A —— B — - Bk 2 + - Bk 2 In kx 4 —— Ak 2 x\ H—— Bk 2 x 2 
2 16 8 4 16 4 128 4 


2x1 


+ 2 /o 


7 11 




120 


47 _ 1 

-B — 3^4x4 —— Bx4 — - Bk 2 Xi In /ex4 + 67Lxl 


1084 

T05~ 


211 

960 


Bxl — Bk 2 xj + Bk 2 xi In kx 4 + 0(xf) + O(y 0 ). 


(92) 


Examining this, we see that to zeroth order in 7/0 the result is in exact agreement with 
© and (|SH) . At y$ order, however, extra terms appear that are not present in (|90|) . 
Nonetheless, at sufficiently small times the effective theory is still valid as these ‘extra’ 
terms are subleading in X 4 : in this limit we find 

TDE 1 

<h 4 =- y + Af + - B 4 k 2 In kx 4 + 0(x 4 ) + 0(/e 4 ) (93) 

X 4 L 

(the superscript E indicating early times), in accordance with the four-dimensional effective 
theory, where 

A f = a ° 4 + ^Byl Bf = S». (94) 

At late times such that x 4 1 (but still on sufficiently long wavelengths that kx 4 -C 1), 
we fold on the positive-tension brane 

4,4 = + + + 

J A B Ak 2 173Bk 2 M: 2 In k Bk 2 lnx 4 

i 3x1 24x1 8x1 960x1 18 x1 12x1 

— - B + - Bk 2 + -Bfc 2 I 11 X 4 + —Bk 2 x^j + O ^^ 3 ^ + O(yo)- (95) 

To zeroth order in y 0 , the results again coincide with the effective theory prediction © 
and However, at y$ order extra terms not present in the four-dimensional effective 

description once more appear. In spite of this, at sufficiently late times the effective descrip¬ 
tion still holds as these ‘extra’ terms are suppressed by inverse powers of X 4 relative to the 
leading terms, which are 

<f >4 = A 4 — —\— ~ A 4 k 2 x\ + Oik 2 In kx 4 ) (96) 

xi 8 
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(where the superscript L indicates late times), in agreement with the four-dimensional ef¬ 
fective theory. (Since X 4 1, we find k <C kx 4 -C 1, and we have chosen to retain terms of 
0(k 2 x \) but to drop terms of 0(k 2 ). The term of 0(x f 2 ) is much larger than 0(k 2 ) and so 
is similarly retained). Fitting this to (1951) . we find j43| 


/\L _ AO _ 
s±4 - ^T-4 


Byo, 


B L 4 = Bl + 



(97) 


Comparing the amplitudes of the two four-dimensional modes at early times, Af and 
Bf, with their counterparts and B% at late times, we see clearly that the amplitudes 
differ at y q order. Using (ITTT1) . we hnd: 


AL A = ( 1 -^o 2 

B L a) \&t l + (| + S)z/o 2 / \Bf) 

Hence the four-dimensional perturbation modes (as defined at very early or very late times) 
undergo mixing. 



VIII. CONCLUSION 

In this paper we have developed a set of powerful analytical methods which, we believe, 
render braneworld cosmological perturbation theory solvable. 

Considering the simplest possible cosmological scenario, consisting of slowly-moving, flat, 
empty branes emerging from a collision, we have found a striking example of how the four¬ 
dimensional effective theory breaks down at first nontrivial order in the brane speed. As 
the branes separate, a qualitative change in the nature of the low energy modes occurs; 
from being nearly uniform across the extra dimension when the brane separation is small, 
to being exponentially localized on the positive-tension brane when the branes are widely 
separated. If the branes separate at finite speed, the localization process fails to keep up 
with the brane separation and the low energy modes do not evolve adiabatically. Instead, a 
given Kaluza-Klein zero mode at early times will generically evolve into a mixture of both 
brane-localized zero modes and excited modes in the late-time theory. From the perspective 
of the four-dimensional theory, this is manifested in the mixing of the four-dimensional 
effective perturbation modes between early and late times, as we have calculated explicitly. 
Such a mixing would be impossible were a local four-dimensional effective theory to remain 
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valid throughout cosmic history: mode-mixing is literally a signature of higher-dimensional 
physics, writ large across the sky. 

As we show in a companion paper 3], this breakdown in the four-dimensional effective 
description has further ramifications for cosmology. A key quantity of interest is the comov¬ 
ing curvature perturbation ( on the positive-tension brane, which is both gauge-invariant 
and, in the absence of additional bulk stresses, conserved on long wavelengths. We show 
that, at first nontrivial order in the brane speed, ( differs from its four-dimensional effec¬ 
tive theory analogue, £ 4 . Hence, while the five-dimensional ( is exactly conserved on long 
wavelengths in the absence of bulk stresses, the four-dimensional effective theory £4 is not 
precisely conserved, contrary to the predictions of the four-dimensional effective theory. This 
has important implications for the propagation of perturbations before and after the bounce 
in cosmologies undergoing a big crunch/big bang transition, such as the ekpyrotic and cyclic 
models. 

The methods developed in this paper should moreover readily extend to other braneworld 
models; for example those containing matter on the branes, and to models deriving from 
heterotic M theory which are better motivated from a fundamental perspective |22j]. A 
'urther application would be to probe the dynamical evolution of braneworld black holes 


and black strings 3 ] in an expanding cosmological background. 

In conclusion, the strength of the expansion about the scaling solution developed in this 
paper lies in its ability to interpolate between very early and very late time behaviors, 
spanning the gap in which the effective theory fails. Not only can we solve for the full 
five-dimensional background and perturbations of a colliding braneworld, but our solution 
takes us beyond the four-dimensional effective theory and into the domain of intrinsically 
higher-dimensional physics. 
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Appendix A: FIVE-DIMENSIONAL LONGITUDINAL GAUGE 
1. Gauge invariant variables 

Starting with the background metric in the form (fTjl. the most general scalar metric 
perturbation can be written as 0 

ds 2 = n 2 (—(1 + 2$) dt 2 — 2W dt dy + t 2 (l — 2T) dy 2 ) — 2Vj« dx l dt + 2f 2 Vj/3 dy dx l 

+b 2 ((1 - 2'k) 8ij - 2V i V j x) dx i dx j . (Al) 

Under a gauge transformation x A —> x A + £ A , these variables transform as 

<i> <i, _ ,c' _ £' r A _ fi/IL 

n n 

r-r + ^-h )? + ?- + ?-, 

t n n 

w ^w-£ t + t 2 i\ 

b 2 ■ 

Oi —> ol — £* 4 —-£ s , 
n z 

X^X + T, (A2) 

where dots and primes indicate differentiation with respect to t and y respectively. Since a 
five-vector has three scalar degrees of freedom £*, and = V*£ s , only four of the seven 
functions (<f>, T, W, a, (3, v h, y) are physical. We can therefore construct four gauge-invariant 
variables, which are 

t , ~ ~ n ~ n' 

^inv = - a - a - 13—, 

n n 

~ 1 Tl ~ T\! 

r inv = r + ^' + - d + d — Y (3 —, 

t n n 

W inv = W — a' + t 2 f3, 
b b' ~ 

dkiiv = ^ + - d + — $, (A3) 

b b 

where d = a - §y and (3 = 0 + ^x'- 

In analogy with the four-dimensional case, we then define five-dimensional longitudinal 
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gauge byy = a = /3 = 0, giving 


$inv = $L, 

P _ P 

1 inv - 1 - L? 

^inv = Wl, 
^inv = ^L, 


(A4) 


i.e. the gauge-invariant variables are equal to the values of the metric perturbations in 
longitudinal gauge. This gauge is spatially isotropic in the x l coordinates, although in 
general there will be a non-zero t-y component of the metric. 


2. Position of branes 

In general, the locations of the branes will be different for different choices of gauge. In 
the case where the brane matter has no anisotropic stresses this is easy to establish. Working 
out the Israel matching conditions, we find that /3 on the branes is related to the anisotropic 
part of the brane’s stress-energy. If we consider only perfect fluids, for which the shear 
vanishes, then the Israel matching conditions give fi(y = ± 2 / 0 ) = 0. 

From the gauge transformations above, we can transform into the gauge a = y = 0 using 
only a £ s and a £* transformation. We may then pass to longitudinal gauge (a — (3 — x — 0) 
with the transformation £ y = (3 alone. Since (3 (and hence /3) vanishes on the branes, £ y must 
also vanish leaving the brane trajectories unperturbed. Hence, in longitudinal gauge the 
brane locations remain at their unperturbed values y = ± 2/0 • Transforming to a completely 
arbitrary gauge, we see that in general the brane locations are given by 

y = ±y 0 - (A5) 
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Appendix B: POLYNOMIAL EXPANSION: BACKGROUND 


Using the expansion in Dirichlet/Neumann polynomials presented in Section IV to solve 
for the background geometry, we find 


«o = I-ii!/U^ji(8-9«' 2 )!/S-^0«(136 + 900i 2 + 375i 4 )9“ 


+ —— t (3968 + 354816 t 2 - 3485441 4 - 36015 1 6 ) y% + O(y 0 10 ) 
Vi = -1 + (4 - 2t 2 ) vl - ^ (61 - 20880 1 2 + 19440l 4 ) yf, 


(61 — 20880 1 2 + 194401 4 ) y 4 


277 743 1 2 6771 4 1011 6 

24192 20~~ + _ 6 3 


vt + o(i/o) 


7\4 3t 3 y 0 2 , t 3 (—28 + 33 1 2 ) Dq , t 3 (1984 - 6776t 2 + 27151 4 ) y 0 6 , /OQ , 

yv 4 = —g—i-^-1-gg- \-u\y o) ( Bd ) 

AT 1 (—7 + 1800 1 2 + 540 t 4 ) yl 

5 “ _ 120 1800 

(323 — 990528 1 2 + 2207520 t 4 + 362880 f 6 ) y$ 6 

201600 + ^ ^ ' 
N e = t 3 y 2 + ^f 3 (-142 + 3711 2 ) y 4 + O(y 6 0 ) (B5) 


+ O(y 0 6 ) 


N 7 = — 


(—9 + 203841 2 + 23520 t 4 ) y 2 


N* = 


5040 
3 t 3 y 2 


94080 


+ O(y 4 0 ) 


+ 0( V q) 


Ng = -+ 0(y o 2 ) 

362880 K 0J 
N w = O(y 2 0 ) 


id 

2 V 8 


-171 2 171 4 55 1 ( 


311 2 9 t 4 2331 6 245t 8 \ „ in , 

W-- + -9CT-W Uo + 0 «) 


-2 t 3 y 2 „ + 


(B10) 


4 yl f—5t 4 9t 6 \ 4 . g 

-f+(—+ —jrf + ofeS) 

, , /7371 3 581 5 \ , . 

Vo+ UnS-5- !/0 + O(!/ 0 6 ) 


Vt + (-^ + - 271 7 ^ yt + 0(yl) (Bll) 


t 4 V 2 

-f + 0(yt) 


(B12) 

(B13) 

(B14) 
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+ 0{yt) 


(B15) 

(B16) 

(B17) 


-131 3 ?/o 2 
60 
0{yl ) 

0{yl ), 

where we have set L = 1 for clarity. (To restore L simply replace t —» f/L). The solution 
has been checked to satisfy all the remaining Einstein equations explicitly. 

Appendix C: POLYNOMIAL EXPANSION: PERTURBATIONS 
1. All wavelengths 

Throughout this Appendix, we shall set the AdS radius L to unity. It is then simple to 
restore L by dimensions where needed, i.e. by setting t —> t/L and k —> kL. Note that the 
coordinate y is dimensionless. 

Using the Dirichlet/Neumann polynomial expansion to solve for the perturbations, the 
solution may be expressed in terms of the original longitudinal gauge variables as 

= v£\y,t)F(°\t) +v£\y,t)FU(t) (Cl) 

where 

F (n \t) = AJ n (kt ) + BY n (kt) (C2) 

for n — 0, 1 and 7 = 0.577... is the Euler-Mascheroni constant. The constants A and B are 
arbitrary functions of k. In order to be consistent with the series expansion in t presented 
in Section III, we must set 

A = 12 A + 2 B k 2 (In 2 - 7 ) - IkM + 233 ^ Vo + 0 (y 6 ) (C 3 ) 

B = Bk 2 n + O(y e 0 ). (C4) 


97 = 
Qs = 
Q9 = 
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The polynomials V^ are then given (for all k and A) by 


-p(o)/. x _ _1 ty_ A 2 (-2y 2 + yg) ty ( lly 2 + 3 (-11 + 3A 2 ) yg) 

* ^ ,V) 6 + 3 + 12 36 

A 2 (-525 y 4 + 90 (19 - 5A 2 ) y 2 y 2 + (511 - 180A 2 + 45 A; 2 A 4 ) yg) 


2160 


A?/ 

2160 


3 ( - 92 + (9 + 4 A; 2 ) A 2 ) y 4 + 30 (92 - (219 + 4 k 2 ) A 2 + A: 2 A 4 ) y 2 y 2 + 


6900 + (20087 + 300 k 2 ) A 2 - 90 (13 + k 2 ) A 4 + 90 A: 2 A 6 ) yg + O(yg), 


(C5) 




1 

A; A 


1 _ Ay A 2 (3 y 2 + (-3 + k 2 A 2 ) y 2 ) 

2 2 + 12 


ty ((3 + At A 2 ) y" + 3 (—3 + (3 — Ar) A 2 + 2 A; 2 A 4 ) yg) 


2160 


36 


+ 


75 (6 + k 2 A 2 ) y 4 + 90 (-12 + 3 (2 - A: 2 ) A 2 + 2 A; 2 A 4 ) y 2 yg + 


(718 + A; 2 A 2 (-101 + 225 A 2 )) yg - 


ty 

2160 


3 (3 + 2 (-9 + 31 k 2 ) A 2 + 2 A; 2 A 4 ) y 4 + 


30 (-3 + (219 - 62 k 2 ) A 2 + 16 A: 2 A 4 ) y 2 yg + 

(225 + (-20104 + 4650 A; 2 ) A 2 + (1215 - 1822 k 2 ) A 4 + 765 k 2 A 6 ) yg 

+O(y 0 6 )- 


(C6) 


Since the F^ are of zeroth order in y 0 , the solution for 4>£ 1° a given order less than O(yg) 
is found simply by truncating the polynomials above. (Should they be needed, results up to 
0(yg 4 ) can in addition be found at j30(). 

In a similar fashion we may express the solution for as 


*l = V^(y,t)F^(t)+vS\y,t)F^(t), 


(C7) 
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where F^> is defined as above and 

wo)/, s _ l_ t _y, g. ( V 2 + Vo) , Ay (~2y 2 + 3 (2 — 3A 2 ) yg) 

* ^ ,?/J 6 3 + 6 36 

A 2 (120 y 4 + 450 (—2 + A 2 ) y 2 yg + (644 + 450A 2 — 45 A; 2 A 4 ) yg) 

2160 + 

^f^3 (—2 + (9 + 4 A; 2 ) A 2 ) y 4 + 30 (2 + 4 (9 - A: 2 ) A 2 + k 2 1 4 ) y 2 yg + 

(-150 + (-2863 + 300 k 2 ) t 2 - 90 (13 + k 2 ) t 4 + 90 k 2 A 6 ) ygj + 0(y°), 

(C8) 


^S\t,y) 


l 

kt 


ty t 2 ( 3 y 2 + (3 + A; 2 A 2 )yg) 


+ 


2 12 
ty ((3 + A; 2 A 2 ) y 2 + 3 (—3 + (3 — k 2 ) t 2 + 2k 2 A 4 ) yg) 


A 2 


2160 


36 


15 (12 + 5 A; 2 A 2 ) y 4 + 90 (-6 - 3 (-2 + k 2 ) A 2 + 2 A; 2 A 4 ) y 2 yg + 


-752 + (540 - 101 k 2 ) A 2 + 360 A; 2 A 4 ) yg + 


^o((9 + 6(-9 + F) + + 

30 (-3 - (33 + 2A; 2 ) A 2 + 7A; 2 A 4 ) y 2 yg + 

(225 + 2 (1288 + 75 k 2 ) A 2 + (1215 - 1012 k' 2 ) A 4 + 765 k 2 A 6 ) yg j 
+O(y 0 6 ). (C9) 


Finally, writing 

W L = V${y,t)F^\t) +p£ ) (y,A)F«(A), (CIO) 

we find 

Purify) = — ^ A 2 ( a / 2 2 / o ) ( — 30 + 30 A y — 25 (y 2 + (-5 + 3 A 2 ) yg) + 

ty (21 y 2 +(-149 + 75 A 2 ) y 2 ) ) +O(y 0 6 ), (Cll) 


"Pwi^v) = ~Q^k t2 ^y 2 ~y^ ( 3 0y-5A; 2 A (2y 2 + (-10+ 3 A 2 ) yg) + 

y ((12 + 11 A; 2 A 2 ) y 2 + (—38 + (60 — 69 A; 2 ) A 2 + 15 A; 2 A 4 ) yg) j + 0(yg). 

(C12) 
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2. Long wavelengths 


On long wavelengths, F^ reduces to 

F m (t) = 12A-?^ + ?^F^ + 0(k 2 )+0( V !), (C13) 

F (1> (t) - (e A t - 0 + - HHlid ~) b) k + 0(fc 2 ) + O(y 0 6 ). (C14) 

For convenience, we provide here a separate listing of the metric perturbations truncated at 

0(k 2 ): 

= (A-^) + (?j- + Aty) + ^ (B(y 2 - Ay 2 ) - AAt 2 (yl + y 2 )) 

+^(5(3y 2 (-4 + t 2 ) + Ay 2 ) + AAt 2 (3y 2 (-19 + 3t 2 ) + 19 y 2 ) 

+ (i/Lt 2 (y 0 4 (29 - 6 1 2 ) + 3y 2 (13 - 2 t 2 )y 2 - 10y 4 ) 

° 4(56 - 45t2) + 1% ° 2( “ 16 + + 100?/4)) ) 

+ ^(2At 2 (5y^(905 - 1338t 2 + 75t 4 ) + 10yg(-181 + 219f 2 )y 2 + 181y 4 ) 

+5(y 4 (50 - 3509t 2 + 135t 4 ) + 5y^(-4 + 235 t 2 )y 2 + (2 - 12f 2 )y 4 )) 

+O(y 0 6 ), (C15) 

A! L = 2A- V -{B + At 2 ) + \{2At 2 (yl + y 2 ) + 5(-y 2 + 2y 2 )) 

-^(4At 2 (3y 2 (-l + 3f 2 ) + y 2 ) + 5(3y 2 (-4 + t 2 ) + 4y 2 )) 

~^48 (^4^ 2 (2yo(17 + 3t 2 ) + 3y^(—7 + 2t 2 )y~ + y 4 ) 

+5(y 4 (8 + 15f 2 ) + 3yo(-8 + 5t 2 )y 2 + 8 y 4 )) 

(27lt 2 (25y 0 4 (l + 42f 2 + 15f 4 ) - 10y 2 (l + 39 t 2 )y 2 + y 4 ) 

+5(y 4 (50 + 721t 2 + 135f 4 ) - 5y^(4 + 47t 2 )y 2 + (2 - 12t 2 )y 4 )) 

+O(y 0 6 ), (C16) 

W L = 6At 2 (-y 2 +y 2 )-t(B + 3At 2 )y(-y 2 + y 2 ) 

+ -0(—yg + y 2 )(—9h>y5 + 20A(y5(—5 + 3t 2 ) + y 2 )) 

-j^hy(-y 2 + y 2 )(120/lt 2 (y 2 (-26 + 9t 2 ) + 3y 2 ) 

+5(y 2 (-152 + 105f 2 ) + 48y 2 )) + O(y 0 6 ). (C17) 

Terms np to 0(yg 4 ) are available at 0 
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Appendix D: PERTURBATIONS FROM EXPANSION ABOUT THE SCALING 
SOLUTION 


Following the method presented in Section V, the perturbations were computed to 0(yf). 
On long wavelengths, the five-dimensional longitudinal gauge variables take the form 

= fo + Vo(fi + ft In (1 + x 4 ) + ft In (1 - x 4 ) + ff In (1 - c jx 4 ), (Dl) 

= ft + Voifi + ft ln (! + % 4 ) + ft In (1 - x 4 ) + ft In (1 - ux 4 ), (D2) 

W L = e~* x * (ff + y 2 0 (ff + ff In (1 + x 4 ) + ff In (1 - x 4 ) + ff ln (1 - ux 4 )) , (D3) 


where the / are rational functions of x 4 and u. For <f+, we have 


fo 


J 1 


ft 


-^-- (16 B - 1§Bujx 4 - 2 ( 8 A + B - ABu 2 )xl + 2(-8 A + 3 B)ux\ 

16^2 (—1 + x 4 ) V 

+(8A-3B)(3 + u; 2 )x^, (D4) 

-f (8Axl ( — 580x 4 + 95^4 — 576u 5 xi + 28 IX 4 + 96o; 6 X4 — 60x1 

960x^(-l +xlf V 4V 

+5o/x 4 (40 + 167x4 ~ 39x 4 ) + 20ca 3 (-19 - 5x 2 + 48x 4 ) 

+10ca 2 x 4 (-78 - 117x 4 + xj + 2x^) + 4ca(285 + 100x^ - 25x^ - 29x^ + 5xf)) 
+B(- 1920 + x 4 (480cu + 80(91 + 15ca 2 )x 4 

-160ca(5 + 4ca 2 )x 2 + 40(-273 - U 6 ca 2 + 7u 4 )xl - 20ca(-649 + 127o; 2 )x^ 
+4(1231 - 3285cu 2 + 2060ta 4 )x| - 4ca(2036 - 2U5ca 2 + 1152u 4 )x% 

+ (1107 + 8102ca 2 - 4905ca 4 + 768o/)xJ + 12ca(233 + 48o/(-5 + 3o/))x® 

-(3131 + 1582a; 2 - 585c/ + 288c/)x^ - 772cjxI° + 20(85 + 29u 2 )xt 
+180a;xf - 120(3 + c/)xf))), (D5) 

--—-- —(b(48 + (36 — 48ca)x 4 + 12(—7 — 6 a; + 2c/)x 2 

48x 4(—1 + x 2 ) ^ 

+4(—13 + 6 a; + 9a/)x 3 + (60 + 80a; — 12c/)x 4 — 4 (—8 — 6 a; + 9u 2 )x^ 

—3(8 + 7a; + 4a/)x 4 + (—11 + 5cu 2 )xl + 3uxf) + 8Axfx 4 + a/x 4 
—a;(l + x 2 )) j, (D 6 ) 
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/.f = --—--^(5(-48 + 12(3 + 4a;)x4-12(-7 + 6co + 2a; 2 )x;| 

48x|(—1 + x|) 

+4(—13 — 6a; + 9a; 2 )x 3 + 4(—15 + 20a; + 3a; 2 )x4 — 4(—8 + 6a; + 9a; 2 )x| 

+3(8 — 7a; + 4a; 2 )x® + (—11 + 5a; 2 )xJ + 3ux\) 

-{-RAx^X/i + a; 2 X 4 — a;(l + x 2 ))^, 

0 _ 35(—1 + UX 4) 2 
h ~ 2 xj(-l + x\f 

Similarly, for we find 


(D7) 

(D8) 


/( 




ft 




/" 3 > 

J 2 

4.41 

■hi 


f 




---7- (165a; — 4(871 + 5(—1 + 2o; 2 ))x 4 + 2(8 A — 3B)lux 4 

16x 4 (—1 + xf) V 

+ (—8A + 35)(—3 + a+)x 3 ^, (D9) 

---—^ ( - 480 Bu> - 2405(—7 + 5co 2 )x 4 + 1605cu(5 + 4ca 2 )^ 

960x 3 (—1 + x|) ' 

-405(143 - 104a; 2 + a o A )x\ + 20co(5(197 - 155ca 2 ) + 87l(-3 + a; 2 ))^ 

+20(-87l(34 - 21u 2 + ca 4 ) + 5(98 - 369w 2 + 101a; 4 ))x 3 
—4a;(2007l(—14 + 5a; 2 ) + 5(34 - 975ca 2 + 288c/))x^ + (4071(55 - 234co 2 + 67co 4 ) 
+5(2455 + 838a; 2 - 1005a/ + 192a/))xJ - 4/35(53 + 120a/ - 36a/) 

+871(155 - 120a; 2 + 36c/))x| + (5(-3515 + 1042a; 2 + 225a/ - 72a/) 

+871(89 + 170a; 2 - 75a/ + 24w 6 ))xl + 4(23271 + 1935)a;xf - 20(871(1 + a/) 


+5(—91 + 29a/))x 41 - 20(871 + 9B)cux\ 2 + 1205(-3 + a/)x 43 ), (DIO) 

-7 (8Axl(x 4 + a/x 4 — a;(l + x 2 A) + 5(36 + 60x 4 — 36x4 — 84x? + 

48x|(-l + xl) 3 V 

24x| + 5x® + u 2 x^(2A + 36x4 — 12x 2 — 36x 3 — 12x4 + 5x 3 ) 

+a;(—48 - 72x 4 + 24x 2 + 80x 3 + 24x 4 - 21x 3 + 3xJ))), (Dll) 


- 7 ( 8-4x5 (—X 4 — a; 2 X 4 + a;(l + x 4 )) + 5(—36 + 6 OX 4 + 36x5 — 

48x 4 (-1 + x 2 ) 3 V 

84x 3 + 24x| - 5x® - a/x 4 (—24 + 36x 4 + 12x 2 - 36x 3 4 + 12x 4 + 5xf) 

+a;(—48 + 72x4 + 24x^ — 80x 3 + 24x^ + 21x| — 3xJ)) j, 

—35(—1 + CUX 4) 2 
2 x |(—1 + x |) 2 


(D12) 

(D13) 
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Finally, for Wl, we have 


w 


fo 

fW 

J 1 


fW _ 

J 2 — 


fW _ 
J 3 — 


fW _ 
J 4 — 


(—1 + u; 2 )a; 4 (—24+4a; 4 (—2 + u;a; 4 ) + B{— 18a; 4 + u;(—8 + 9a; 2 ))) 

8(—1 + a;|) 2 

48 o 1 k - + 1+^( 8 ^ (1500 + 84 ^ 1( " 12 + ^ 

—6o; 2 (50 + 100.x- 2 - 427x1 + x\) + 3ua 3 a; 4 (60 + 585^ - 160x1 + 3x1) 
+co 5 (168a;| - 36^) - 6^(-590 + 243a^ + 201x^ - Ax 6 4 + 10a?|) 
+coa; 4 (—1560 - 4935a; 2 + 2000x^ - 265a:® + 92a;®)) 

+B (1440 + a; 4 (—84u; 4 a;|(24 + 28a; 2 + 3x1) 

+12u 5 x 6 4 (4:0 + 6a;^ + 9x1) + 6oa 2 a;^(-450 + 964a;^ + 863a;^ + 3x1) 

—3a; 3 x 2 (—40 - 748a; 2 - 2333a; 4 + 672a;® + 9a;®) 

—6a; 4 (2160 - 5490a;^ + 3770a^ - 2249a;^ + 597x^ - 588a;^ 0 + 90x\ 2 ) 
+o;(—2160 + 18920a; 2 - 53216a4 + 20629a;^ - 11216x1 + 5579a;^° 


14 


-2236a;4 + 360a; 4 


( 1 -^) 
24(a; 4 — a; 3 ) 4 


(b( - 144 + 10S(—1 + 2u)x 4 - 24(-5++)(3 + 2+)+ 


(D14) 


(D15) 


—36(—11 + +(16 + u))x\ + 24+(-29 + 7+)+J 

+36(—2 + u (—2 + 7co))xl + 3(1 + u;)(3 + 32u;)a: 4 + 7u;(l + uj)x 7 4 

+27(1 + uj)x\ — 9ua(l + u>)x^j + 24+4(1 + ui)a;®(—1 — 3a; 2 + co(a; 4 + a; 3 ))^ 

- ^ - JB( - 144 + 108(1 + 2ca)a; 4 - 24(5 + cu)(-3 + 2io)a; 2 

24x 4 (—1 + x 4 ) ' 

+36(—11 + (-16 + o;)o;)a; 4 + 24oa(29 + 7u)x 4 4 - 36(-2 + u{2 + 7cu))x 5 4 
+3(—1 + u >)(—3 + 32u;)a;4 — 7(—1 + u)ux 7 4 — 27(—1 + a;)x® 

+9(—1 + a;)u;xl) — 24+4(—1 + u)x 4 (—1 — 3a; 4 + cu(a; 4 + 

3B(—1 + u;a: 4 ) 2 (4 — 10a;| + cua; 4 (—1 + 7x 4 )) 


,(D16) 


x: 


x 4 (—1 + x 2 4 ) 4 

Results including the k 2 corrections can be found at [301- 


(017) 

(018) 
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Appendix E: BULK GEODESICS 


To calculate the affine distance between the branes along a spacelike geodesic we must 
solve the geodesic equations in the bulk. Let us first consider the situation in Birkhoff-frame 
coordinates for which the bulk metric is static and the branes are moving. The Birkhoff- 
frame metric takes the form 132] 


ds 2 = dK 2 - N 2 (Y) dT 2 + A 2 (Y) d£ 2 , 


(El) 


where for Schwarzschild-AdS with a horizon at Y = 0, 

A 2 (Y cosh(2F/L) 2 cosh(21'o/L) / sinh(2F/L) \ 2 

1 } cosh(2 Y Q /Ly [ ’ cosh(2 Y/L) \smh(2Y 0 /L)J 


(E2) 


At T = 0, the Y -coordinate of the branes is represented by the parameter Y 0 : their subse¬ 
quent trajectories Y±{T) can then be determined by integrating the Israel matching condi¬ 
tions, which read tanh (2 Y±/L) = ±y/l — V± , where V± = (dY±/dT)/N(Y±) are the proper 
speeds of the positive- and negative-tension branes respectively. From this, it further follows 
that Y 0 is related to the rapidity y 0 of the collision by tanh y 0 = sech(21j)/L). 

For the purpose of measuring the distance between the branes, a natural choice is to use 
spacelike geodesics that are orthogonal to the four translational Killing vectors of the static 
bulk, corresponding to shifts in x and T. Taking the x and T coordinates to be fixed along 
the geodesic then, we find that Y]a is constant for an affine parameter A along the geodesic. 

To make the connection to our original brane-static coordinate system, recall that the 
metric function b 2 (t,y ) = A 2 (F), and thus 


V 2 _ (+ bb iV y, A ) 2 
’ A b 4 - e 2 


= n 2 (—t\ + t 2 y 2 x ), 


(E3) 


where we have introduced the constant 6 = tanh y 0 = V/c. Adopting y now as the affine 
parameter, we have 


0 = {b 2 t b 2 + n 2 (6 4 - 9 2 ))t 2 y + 2b f b, y b 2 t, y + ( b 2 y b 2 - n 2 t 2 (6 4 - 6 2 )), (E4) 

where t is to be regarded now as a function of y. 

We can solve this equation order by order in y 0 using the series ansatz 

OO 

t(y) = ^c n y n , (E5) 

71=0 
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where the constants c n are themselves series in t/q- Using the series solution for the back¬ 
ground geometry given in Appendix B, and imposing the boundary condition that t(y 0 ) = t 0 , 
we obtain 


c 0 


Cl 

c 2 

c 3 

c 4 

C 5 

C 6 

C 7 


to + ^ - 2 tlvl + (t ° + - tl (1 + 5«g) vl + f At + m ° 


2 ~ Ut7U ' 24 

tg (13 + 250tg + 795tg) i/ 0 7 


720 


+ 4«8 y§ 


60 


+ 0(4) 


2t^„ 2 +(UI + 5i p y‘-8« + 


9U| 

180 


+ 


t 0 t 0 (l + 6tg)4 


+ *o Vo ~ ~ 2i o + 4 *o ) Vo + 


2 4 

5tg4 *o (25 + 201 tg) 4 


23 tp 

~6~ 

5 \ „ ,4 


+ 


53 tg 


4 + 0(4) 


(E6) 

(E7) 


51 0 


+ 


24 

61 *o4 

60 
611 0 

720 
0 + O(y 0 ). 


18 

5t ° , l^o\ 2 

48 + 4 ) Vo 


48 


+ 0(4) 


(to + 23 tg) 1 /g 


5*0 4 

12 


+ 0(4) 


+ 0(4) 
+ 0(4) 


+ 0(4) (E8) 

(E9) 

(E10) 

(Ell) 

(E12) 

(E13) 


Substituting t 0 = x 0 /|/o and y = cm/ 0 we End a;(a;) = ^o/t/o + O (y 0 ), i.e. to lowest order in y 0 , 
the geodesics are trajectories of constant time lying solely along the u direction. Hence in 
this limit, the affine and metric separation of the branes (defined in (1HU1B must necessarily 
agree. To check this, the affine distance between the branes is given by 


da 

T 


ry o 


n 


' -y o 


Vt 2 — t ' 2 d y 

(t 0 + 5 tg) 4 , 4 , (*o - 10 tg + 159 tg) yl 


~ 4 1 0 y 0 + 


2toyo+ 3 60 

(to + 31115 tg - 5523 tg + 12795 tg) y 7 0 


2520 


+ 0(4), 


(E14) 

2 (tp + 30 tg) yt 
3 

(E15) 


which to lowest order in y 0 reduces to 

si ^ ^ zv»5 

T —2*0 + — + — 


+ 


853 4 
168 


+ 0 ( 4 ) + 0 (?/ 5 ), 


(E16) 


in agreement with the series expansion of (pH) . (Note however that the two distance measures 
differ nontrivially at order 4)- 
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To evaluate the perturbation 8d a in the affine distance between the branes, consider 


8 / sJg pu xPx v dA = 


dA 


^/g pr7 xPx c 


{8g pu x^x v + g ll „ )K 5x K x ll x'' + 2 g pv x^8x v 


x,,8x v 


sJg pa xPx° 


1 r iX 

2 J v /g pa xPx a 


(E17) 


where dots indicate differentiation with respect to the affine parameter A, and in going to the 
second line we have integrated by parts and made use of the background geodesic equation 
x a = ^g p ,u,aX tl x l/ and the constraint = 1. If the endpoints of the geodesics on the 

branes are unperturbed, this expression is further simplified by the vanishing of the surface 
term. Converting to coordinates where t 0 = Xo/yo and y = uyo, to lowest order in y 0 the 
unperturbed geodesics lie purely in the u direction, and so the perturbed affine distance is 
once again identical to the perturbed metric distance EH)- 
Explicitly, we find 


8da = _ 2 (B + Atl ) y 0 _ f B (4 + 3t 1 2 0 ) A (t 0 + 9t 3 0 ) \ 3 + (_ AB + AAt 2) 4 
L t 0 V 12t o 3 J y ° V 0)y ° 

B (2+ 2169+)+ 135+)) +2 At 2 0 (1 + 1110tg + 375^ 5 

120+) ) y ° 

4 At q (1 + 42^5) — B (4 + 57fg) A 6 
6 ) Vo 

B (4 + 88885+) + 952866+) + 28875+)) + 4 At 2 0 (1 - 152481+) + 293517+) + 36015+)) 


Vo 


10080+ 


+0(% 8 ), 


(E18) 


which, substituting t 0 = Xo/yo and dropping terms of 0(?/g), reduces to 


8 d a 
L 


2 B „ B , , 9 ~ , 25 , , 275 
■—-2Ax 0 -jX 0 -3Ax 0 --Bx 0 -—Ax 0 - — 


Bxl 


343 
1 24 


Axq+0(x 0 ), (E19) 


where B = By$. Once again, this expression is in accordance with the series expansion of 
EH)- However, the perturbed affine and metric distances do not agree at 0(yl). 
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